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ABSTRACT 


This  thesis  investigates  the  simulation  performance  of  multiple-input  multiple- 
output  (MIMO)  systems  utilizing  Alamouti-based  space-time  block  coding  (STBC)  tech¬ 
nique.  The  MIMO  communication  systems  using  STBC  technique  employing  both  sin¬ 
gle-carrier  modulation  and  orthogonal  frequency  division  multiplexing  (OFDM)  are 
simulated  in  Matlab.  The  physical  layer  part  of  the  IEEE  802.16a  standard  is  used  in  con¬ 
structing  the  simulated  OFDM  schemes.  Stanford  University  Interim  (SUI)  channel  mod¬ 
els  are  selected  for  the  wireless  channel  in  the  simulation  process.  The  performance  re¬ 
sults  of  the  simulated  MIMO  systems  are  compared  to  those  of  conventional  single  an¬ 
tenna  systems. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  OBJECTIVE  AND  METHODOLOGY . 2 

B.  RELATED  RESEARCH . 2 

C.  ORGANIZATION  OF  THE  THESIS . 2 

II.  MIMO  SYSTEMS  AND  ALAMOUTI  SCHEME . 5 

A.  MULTIPLE-INPUT  MULTIPLE- OUTPUT  (MIMO)  SYSTEMS . 5 

1.  Single-Input  Single-Output  (SISO)  Systems . 5 

2.  Single-Input  Multiple-Output  (SIMO)  Systems . 6 

3.  Multiple-Input  Single-Output  (MISO)  Systems . 7 

4.  Multiple-Input  Multiple- Output  (MIMO)  Systems . 7 

B.  ALAMOUTI  SCHEME . 8 

1.  Maximal  Ratio  Combining  (MRC)  Scheme  for  1x2  Systems . 9 

2.  Alamouti  Scheme  for  2x1  Systems . 10 

3.  Alamouti  Scheme  for  2x2  Systems . 11 

C.  SUMMARY . 13 

III.  MIMO  CHANNEL  MODELS . 15 

A.  STANFORD  UNIVERSITY  INTERIM  (SUI)  CHANNELS  FOR  THE 

IEEE  802.16a  STANDARD . 15 

B.  GENERATING  CORRELATED  MIMO  CHANNELS . 17 

1.  Uncorrelated  MIMO  Channels . 17 

2.  Correlated  MIMO  Channels . 18 

3.  LOS  Component . 19 

C.  GENERATING  CORRELATED  FREQUENCY-SELECTIVE 

MIMO  CHANNELS . 20 

D.  SUMMARY . 21 

IV.  ALAMOUTI-BASED  SCHEMES  OVER  FREQUENCY-SELECTIVE 

CHANNELS . 23 

A.  SPACE-TIME  BLOCK  CODING-SINGLE  CARRIER  (STBC-SC) 

SYSTEMS  OVER  FREQUENCY-SELECTIVE  CHANNELS . 24 

1.  SISO  Systems  over  Frequency- Selective  Channels . 24 

2.  SIMO  Systems  over  Frequency- Selective  Channels . 24 

3.  MISO  Systems  over  Frequency-Selective  Channels . 26 

4.  MIMO  Systems  over  Frequency- Selective  Channels . 27 

B.  SPACE-TIME  BLOCK  CODING-ORTHOGONAL  FREQUENCY 
DIVISION  MULTIPLEXING  (STBC-OFDM)  SYSTEMS  OVER 

FREQUENCY-SELECTIVE  CHANNELS . 29 

1.  SISO-OFDM  Systems  over  Frequency-Selective  Channels . 29 

a.  Channel  Encoding . 30 

b.  Symbol  Mapping . 30 

c.  Inverse  Fast  Fourier  Transform  (IFFT)  Operation . 31 

vii 


d.  Cyclic  Prefix  Addition . 33 

e.  Digital-to- Analog  Conversion  and  RF  Modulation . 34 

f  RF  Demodulation  and  Digital  to  Analog  Conversion . 34 

g.  Cyclic  Prefix  Removal  and  FFT  Operation . 34 

h.  Demapping  and  Decoding . 35 

2.  SIMO-OFDM  Systems  over  F requency-Selective  Channels . 35 

3.  MISO-OFDM  Systems  over  F requency-Selective  Channels . 36 

4.  MIMO-OFDM  Systems  over  Frequency- Selective  Channels . 38 

C.  SUMMARY . 39 

Y.  SIMULATION  RESULTS . 41 

A.  PERFORMANCE  OF  SYSTEMS  USING  ALAMOUTI  SCHEMES . 41 

B.  PERFORMANCE  OF  MIMO-SC  SYSTEMS  OVER  SUI 

CHANNELS . 43 

C.  PERFORMANCE  OF  MIMO-OFDM  SYSTEMS  OVER  SUI 

CHANNELS . 49 

D.  SUMMARY . 54 

VI.  CONCLUSION . 57 

A.  SUMMARY  OF  THE  WORK  DONE . 57 

B.  SIGNIFICANT  RESULTS  AND  CONCLUSIONS . 57 

C.  SUGGESTIONS  FOR  FUTURE  WORK . 58 

APPENDIX  A.  SUI  CHANNEL  PARAMETERS . 61 

APPENDIX  B.  MATLAB  CODE  EXPLANATION . 63 

A.  SIMULATION  PARAMETERS . 63 

B.  CONDUCTED  ITERATIONS . 65 

C.  EXPLANATIONS  OF  THE  SUB-FUNCTIONS . 66 

LIST  OF  REFERENCES . 69 

INITIAL  DISTRIBUTION  LIST . 73 


viii 


LIST  OF  FIGURES 


Figure  1.  Schematic  of  a  SISO  System  (After  Ref.  2.) . 5 

Figure  2.  Schematic  of  a  SIMO  system  (After  Ref.  2.) . 6 

Figure  3.  Schematic  of  a  MISO  System  (After  Ref.  2.) . 7 

Figure  4.  Schematic  of  a  MIMO  System  (After  Ref.  2.) . 8 

Figure  5.  Maximal  Ratio  Combining  for  a  1x2  System  (After  Ref.  4.) . 10 

Figure  6.  Alamouti  Scheme  for  a  2x2  System  (After  Ref.  4.) . 10 

Figure  7.  Alamouti  Scheme  for  a  2x1  System  (After  Ref.  4.) . 12 

Figure  8.  The  Frequency  Response  Plots  of  Different  SUI  Channel  Models . 17 

Figure  9.  Conceptual  Diagram  for  Generating  a  Correlated  MIMO  Channel  Based 

on  SUI  Channel  Models . 22 

Figure  10.  Schematic  of  a  1x2  SIMO-SC  System  Utilizing  MRC  over  Frequency- 

Selective  Channels . 25 

Figure  11.  Schematic  of  a  2x1  STBC-SC  System  over  Frequency-Selective  Channels  ..26 

Figure  12.  Schematic  of  a  2x2  STBC-SC  System  over  Frequency-Selective  Channels  ..28 

Figure  13.  Block  Diagram  of  a  SISO-OFDM  Scheme . 30 

Figure  14.  Gray  coded  (a)  BPSK,  (b)  QPSK  and  (c)  16-QAM  constellations  (From 

Ref.  10.) . 31 

Figure  15.  Allocation  of  Subcarriers  before  the  IFFT  Block  (see  Table  6) . 32 

Figure  16.  OFDM  Symbol  Time  Structure  (From  Ref.  10.) . 33 

Figure  17.  Schematic  of  a  1x2  SIMO-OFDM  System  Utilizing  MRC . 36 

Figure  18.  Schematic  of  a  2x1  MISO-OFDM  System  Utilizing  STBC . 37 

Figure  19.  Schematic  of  a  2x2  MIMO-OFDM  System  Utilizing  STBC . 39 

Figure  20.  The  BER  Performance  Comparison  of  the  Alamouti  Schemes  to  the  MRC 

Scheme  over  Rayleigh  Fading  Channels  (The  simulation  parameters  are 

given  in  Table  7) . 43 

Figure  21.  SUI  Channel  Impulse  Responses  on  a  Sampling  Grid  with  a  Sampling 

Frequency  fs  =  22.857  MHz . 44 

Figure  22.  Schematic  of  Generation  of  192-symbol  Blocks . 45 

Figure  23.  The  BER  Performance  of  the  SISO-SC  System  over  Different  SUI 

Channels  (The  simulation  parameters  are  given  in  Table  8) . 46 

Figure  24.  Schematic  of  Generation  of  192-symbol  Blocks  for  each  Transmit 

Antenna . 47 

Figure  25.  The  BER  Performance  Comparisons  of  Single-Carrier  Schemes  over  the 

SUI-2  Channel  (The  simulation  parameters  are  given  in  Table  9) . 48 

Figure  26.  The  BER  Performance  of  the  MISO  Single-Carrier  System  over 

Correlated  Channels  (The  simulation  parameters  are  given  in  Table  10) . 50 

Figure  27.  The  BER  Comparison  of  SISO-OFDM  with  Various  CP  Sizes  to  SISO-SC 

over  the  SUI-2  Channel  (The  simulation  parameters  are  given  in  Table  12). .53 
Figure  28.  The  BER  Comparison  of  all  Schemes  over  the  SUI-2  Channel  (The 

simulation  parameters  are  given  in  Table  13) . 55 

Figure  29.  A  Conceptual  Diagram  of  the  Iterations  for  the  Outer  Functions . 66 


IX 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


x 


LIST  OF  TABLES 


Table  1.  Transmission  Sequence  of  the  Alamouti  Scheme . 1 1 

Table  2.  Notation  for  the  Received  Signals  in  2x2  Alamouti  Scheme . 12 

Table  3.  SUI  Channel  Characterization  (After  Ref.  13.) . 16 

Table  4.  Transmission  Sequence  for  STBC . 27 

Table  5.  The  Sequence  and  the  Location  of  the  Received  Signals  for  a  2x2  STBC 

System . 29 

Table  6.  Allocation  of  OFDM  Subcarriers  in  the  IEEE  802.16a  Standard . 32 

Table  7.  Parameters  for  Alamouti  and  MRC  Single-Carrier  Simulations . 42 

Table  8.  Parameters  for  SISO  System  Simulations  over  SUI  Channels . 45 

Table  9.  Parameters  for  Simulations  of  the  Single-Carrier  Systems  over  SUI-2 

Channel . 47 

Table  10.  Parameters  for  Simulations  of  the  MISO  System  over  Correlated  SUI-2 

Channels . 49 

Table  11.  OFDM  Symbol  Durations  for  Various  Guard  Ratios . 51 

Table  12.  Parameters  for  Simulations  of  the  SISO-SC  and  SISO-OFDM  Systems 

over  the  SUI-2  Channel  Model . 52 

Table  13.  Parameters  for  Simulations  of  the  all  Single-Carrier  and  OFDM  Systems 

over  the  SUI-2  Channel . 54 

Table  14.  Parameters  of  the  SUI  Channel  Models . 62 

Table  15.  The  Variables  Representing  the  Simulation  Parameters . 64 


xi 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


LIST  OF  ACRONYMS  AND  ABBREVIATIONS 


AWGN 

BER 

CSI 

dB 

eft 

IFFT 

ISI 

MIMO 

MISO 

MRC 

OFDM 

PSK 

QAM 

RF 

SC 

SIMO 

SISO 

STBC 

SUI 


Additive  White  Gaussian  Noise 
Bit  Error  Rate 
Channel  State  Information 
decibel 

Fast  Fourier  Transform 

Inverse  Fast  Fourier  Transform 

Inter  Symbol  Interference 

Multiple  Input  Multiple  Output 

Multiple  Input  Single  Output 

Maximal  Ratio  Combining 

Orthogonal  Frequency  Division  Multiplexing 

Phase  Shift  Keying 

Quadrature  Amplitude  Modulation 

Radio  Frequency 

Single  Carrier 

Single  Input  Multiple  Output 
Single  Input  Single  Output 
Space  Time  Block  Coding 
Stanford  University  Interim 


xiii 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xiv 


ACKNOWLEDGEMENTS 


I  thank  my  wife  Cigdem  for  being  in  my  life  and  for  the  support  and  love  she  has 
given  during  my  studies  at  NPS. 

Thanks  go  out  as  well  to  Ruth  and  Bobby  Moorhatch.  The  support  and  love  you 
gave  us  will  never  be  forgotten. 

Finally,  I  want  to  thank  my  advisor  Prof.  Murali  Tummala  for  his  guidance  and 
patience  during  the  development  of  this  thesis. 


xv 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xvi 


EXECUTIVE  SUMMARY 


Future  military  and  commercial  wireless  communication  systems  are  required  to 
provide  higher  data  rates  and  reliable  communication.  Two  major  challenges  in  system 
design  are  the  limited  spectrum  and  the  fading  caused  by  multipath  components  in  the 
wireless  channel.  To  obtain  higher  data  rates,  larger  bandwidths  are  required.  To  keep  the 
same  reliability  under  severe  channel  conditions,  a  higher  transmitted  power  level  is  nec¬ 
essary.  The  emerging  multiple-input  multiple-output  (MIMO)  communication  technolo¬ 
gies  have  the  potential  to  improve  the  performance  without  increasing  the  bandwidth  or 
the  transmitted  power.  MIMO  systems  exploit  spatial  diversity  by  employing  multiple 
antennas  at  either  side  of  the  communication  link.  MIMO  systems  may  add  robustness  to 
future  military  communication  systems  under  battlefield  conditions. 

Due  to  its  robustness  under  frequency  selective  channels,  orthogonal  frequency 
division  multiplexing  (OFDM)  has  been  adopted  in  several  wireless  communication  stan¬ 
dards,  such  as  the  IEEE  802.11a  local  area  network  (LAN)  standard  and  the  IEEE 
802.16a  metropolitan  area  network  (MAN)  standard. 

The  combination  of  OFDM  and  MIMO  technologies  is  a  potential  candidate  for 
future  wireless  systems.  When  the  spatial  diversity  gain  provided  by  a  MIMO  system  is 
added  to  robustness  to  fading  provided  by  the  OFDM  systems,  the  performance  of  the 
resulting  system  is  significantly  enhanced  compared  to  the  conventional  single  antenna 
and  single-carrier  systems. 

The  main  objective  of  this  thesis  was  to  investigate  MIMO  and  MIMO-OFDM 
systems.  The  space  coding  techniques  and  the  decoding  algorithms  studied  in  this  thesis 
are  based  on  the  widely  accepted  Alamouti  scheme.  The  Alamouti  scheme  was  originally 
proposed  for  flat  fading  channels.  The  extended  scheme,  space-time  block  coding 
(STBC),  has  been  developed  for  frequency- selective,  multipath  channels.  The  STBC 
technique  is  investigated  for  systems  employing  both  single-carrier  and  OFDM  modula¬ 
tion.  Systems  with  two  transmit  antennas  and  a  single  receive  antenna  are  the  most  inves¬ 
tigated  schemes  utilizing  the  Alamouti  and  the  STBC  MIMO  technique.  By  extending  the 

xvii 


existing  results  for  this  case,  in  this  thesis,  we  developed  the  input-output  relations  and 
decoding  equations  for  the  other  possible  antenna  combinations,  such  as  one  transmitting 
and  two  receiving  antennas,  and  two  transmitting  and  two  receiving  antennas. 

The  widely  used  Stanford  University  Interim  (SUI)  channel  models  were  selected 
to  represent  the  wireless  channel.  The  use  of  multiple  antennas  at  both  ends  of  a  commu¬ 
nication  system  results  in  multiple  channels  between  the  transmitter  and  the  receiver.  The 
multiple  channels  of  MIMO  systems  are  constructed  using  the  SUI  channel  models.  In 
practice,  the  spacing  between  antennas  and  the  environment  effects  (angle  of  arrival  and 
angle  of  departure  of  the  electromagnetic  waves)  may  cause  correlation  between  these 
multiple  channels.  We  investigated  the  spatial  correlation  effects  on  both  frequency-flat 
and  frequency- selective  MIMO  channels. 

The  MIMO  systems  with  single-carrier  and  OFDM  modulation  were  simulated  in 
Matlab.  Simulation  results  are  presented  in  the  form  of  the  bit  error  rate  (BER)  curves.  A 
performance  comparison  among  the  various  MIMO  systems  is  reported. 

The  most  significant  result  of  the  simulations  was  that  the  systems  employing 
higher  number  of  antennas  at  either  side  of  the  communication  link  performed  better  than 
the  ones  with  fewer  antennas.  In  the  work  reported  here,  the  best  performing  system  was 
the  MIMO-OFDM  system  with  two  transmitting  and  two  receiving  antennas.  In  a  simula¬ 
tion  using  the  SUI-2  channel  model,  at  a  BER  of  1 0  4 .  MIMO-OFDM  performed  12  dB 
better  than  the  conventional  single  antenna  system  employing  single -carrier  modulation; 
the  total  transmitted  power  was  kept  at  the  same  level  in  each  scheme. 


I.  INTRODUCTION 


One  of  the  biggest  challenges  in  wireless  communication  is  to  operate  in  a  time- 
varying  multipath  fading  environment  under  limited  power  constraints.  The  other  chal¬ 
lenge  is  the  limited  availability  of  the  frequency  spectrum.  Future  commercial  and  mili¬ 
tary  wireless  systems  will  be  required  to  support  higher  data  rates  with  reliable  communi¬ 
cation  under  spectrum  limitations  and  multipath  fading  environments.  Military  communi¬ 
cation  systems  must  maintain  reliable  communication  under  the  conditions  of  hostile 
jamming  and  other  interference  without  increasing  emitted  power  or  requiring  larger 
bandwidth. 

In  order  to  improve  the  reliability  without  increasing  the  emitted  power,  time,  fre¬ 
quency  or  space  diversity  could  be  exploited.  In  time  diversity,  the  received  signal  is 
sampled  at  a  higher  rate,  thus  providing  more  than  one  sample  per  transmitted  symbol.  In 
frequency  diversity,  the  same  information  is  sent  over  a  number  of  carriers  [1].  Both  di¬ 
versity  techniques  require  larger  bandwidth.  To  exploit  space  diversity,  the  same  infor¬ 
mation  is  transmitted  or  received  through  multiple  antennas.  Employing  multiple  anten¬ 
nas  at  the  receiver  and/or  the  transmitter  improves  the  quality  of  a  wireless  communica¬ 
tion  link  without  increasing  the  transmitted  power  or  bandwidth  [2].  Therefore,  the  design 
and  implementation  of  multiple-input  multiple-output  (MIMO)  communication  systems 
is  an  attractive  research  area. 

Orthogonal  frequency  division  multiplexing  (OFDM)  is  a  widely  used  method  in 
wireless  communication  systems.  Due  to  its  effectiveness  in  multipath  channel  condi¬ 
tions,  OFDM  has  been  adopted  by  several  wireless  communication  standards,  such  as  the 
IEEE  802.11a  local  area  network  (LAN)  standard  and  the  IEEE  802.16a  metropolitan 
area  network  (MAN)  standard.  The  combination  of  OFDM  and  MIMO  systems  presents 
better  solutions  by  adding  more  diversity  gain  to  the  conventional  OFDM  systems  em¬ 
ploying  a  single  antenna  at  both  the  receiver  and  the  transmitter  [3].  The  robustness  to 
fading  provided  by  OFDM  is  enhanced  by  the  spatial  diversity  of  MIMO  systems,  and 
the  resulting  performance  of  MIMO-OFDM  systems  is  significantly  improved. 
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A.  OBJECTIVE  AND  METHODOLOGY 

The  main  objective  of  this  thesis  was  to  investigate  MIMO  and  MIMO-OFDM 
systems  and  compare  their  performance  to  the  conventional  single  antenna  systems.  The 
first  step  to  achieve  this  goal  was  to  study  the  fundamentals  of  MIMO  systems  by  inves¬ 
tigating  their  performance  in  a  communication  system.  In  this  thesis,  the  published  tech¬ 
niques  for  the  systems  with  single-carrier  modulation  and  with  OFDM  were  investigated. 
Several  communication  systems  with  various  numbers  of  antennas  utilizing  both  single¬ 
carrier  modulation  and  OFDM  were  developed  in  Matlab.  The  developed  systems  were 
simulated  using  the  widely  used  Stanford  University  Interim  (SUI)  channel  models. 

B.  RELATED  RESEARCH 

Due  to  their  efficiency  in  providing  an  improved  performance  without  increasing 
the  bandwidth  or  the  emitted  power,  MIMO  systems  are  the  subject  of  considerable  re¬ 
search  effort.  Probably  the  most  attractive  scheme  from  the  stand  point  of  implementation 
and  performance  is  the  Alamouti’s  transmit  diversity  scheme  [4]  based  on  maximal  ratio 
combining  (MRC)  [5].  Numerous  studies  have  been  performed  to  investigate  its  perform¬ 
ance  and  led  to  the  development  of  several  variations  [6,  7].  The  scheme  was  originally 
developed  for  flat  fading  channels.  Subsequently,  it  has  been  extended  to  frequency  se¬ 
lective  channel  cases  [8,  9]  and  renamed  as  “the  space-time  block  coding  (STBC)  tech¬ 
nique”.  Later,  the  extended  scheme  has  been  selected  to  be  a  part  of  the  IEEE  802.16a 
standard  [10]. 

Channel  modeling  is  another  research  area  in  MIMO  systems.  The  IEEE  802.16 
Broadband  Wireless  Access  Working  Group  proposed  SUI  channel  models  for  system 
simulation  [11].  The  presented  channel  models  have  been  widely  used  by  many  research 
studies. 

C.  ORGANIZATION  OF  THE  THESIS 

This  thesis  is  organized  into  six  chapters.  Chapter  II  introduces  the  input-output 
relations  in  MIMO  systems  and  Alamouti’s  transmit  diversity  scheme  designed  for  flat 
fading  channels.  Chapter  III  introduces  the  SUI  channel  models.  The  uncorrelated  MIMO 
channels  based  on  SUI  modeling  and  their  extension  to  correlated  MIMO  channels  are 
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also  discussed.  Chapter  IV  introduces  the  STBC  scheme  that  uses  both  single-carrier 
modulation  and  OFDM.  Chapter  V  presents  the  simulation  results  of  the  communication 
systems  developed  in  Matlab.  Chapter  VI  provides  a  summary  of  the  work,  the  conclu¬ 
sions  and  suggestions  for  future  studies. 

Appendix  A  lists  the  SUI  channel  parameters.  Appendix  B  provides  an  explana¬ 
tion  of  the  Matlab  simulation. 
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II.  MIMO  SYSTEMS  AND  ALAMOUTI  SCHEME 


The  goal  of  this  chapter  is  to  introduce  a  widely  used  space  coding  technique 
called  the  Alamouti  scheme.  An  overview  of  the  input-output  relations  of  MIMO  systems 
is  presented  to  describe  the  space  coding  algorithm  at  the  transmitter  and  the  decoding 
process  at  the  receiver  in  the  Alamouti  scheme. 

A.  MULTIPLE-INPUT  MULTIPLE-OUTPUT  (MIMO)  SYSTEMS 

In  order  to  develop  the  input-output  relations  of  single-input  multiple-output 
(SIMO),  multiple-input  single-output  (MISO)  and  MIMO  systems,  the  single-input  sin¬ 
gle-output  (SISO)  system  is  described  first. 

1.  Single-Input  Single-Output  (SISO)  Systems 

The  schematic  diagram  of  a  SISO  system  is  shown  in  Figure  1.  The  time-variant 
channel  impulse  response  from  the  transmitter  to  the  receiver  is  denoted  as  h(r,t)  where 

it  represents  the  response  at  time  t  with  an  impulse  applied  at  time  t  —  T,  which  indicates 
that  in  a  general  setup,  the  system  is  time-varying.  The  input-output  relation  for  a  SISO 
system  is  given  by 

r(t)  =  [  '  h(r,t) x(t -t)cIt  =  h(r,t)* x(t)  (2.1) 

J  0 

where  x(t)  is  the  transmitted  signal  and  r(l)  denotes  the  received  signal  at  time  1 .  The 
upper  limit  of  the  integral  Tlotal  is  the  duration  of  the  impulse  response.  The  operator  * 
denotes  the  convolution  operation  [2] . 


Figure  1.  Schematic  of  a  SISO  System  (After  Ref.  2.) 
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2.  Single-Input  Multiple-Output  (SIMO)  Systems 

A  SIMO  system  has  a  single  transmit  antenna  and  multiple  (Mr )  receive  anten¬ 
nas.  Figure  2  illustrates  the  schematic  of  a  SIMO  system. 


Figure  2.  Schematic  of  a  SIMO  system  (After  Ref.  2.) 

The  channels  ht{ T,t),  i  =  1,2, .... Mr ,  represent  the  impulse  responses  from  the 
single  transmit  antenna  to  the  Mr  receive  antennas.  The  received  signals  at  the  respective 
receive  antennas  are  given  by 

r1(t)  =  hl(T,f)*x(t) 
r2{t)  =  h2{T,t)*x{t) 

rMr  {t)  =  hMr{^t)*X{t)- 

Representing  the  signals  at  the  receiver  antennas  in  a  vector  form,  we  have  the  received 
signal  vector 

r(0  =  [li(0  r2  (0  •••  V(r)]7  -  (2-3) 

Similarly,  the  channel  vector  is  given  by 

1>M=[M»)  K (<)  •••  (2-4) 

Therefore,  we  can  express  the  input-output  relationship  of  a  SIMO  system  as 

r(t)  =  h(r,t)  A.r(t)  (2.5) 

where  the  operator  A  denotes  the  element-by-element  convolution. 
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3.  Multiple-Input  Single-Output  (MISO)  Systems 

Figure  3  depicts  a  simple  MISO  system  with  multiple  ( M t )  transmit  antennas  and 

a  single  receive  antenna.  Input-output  relationship  of  a  MISO  system  can  be  developed 
on  the  lines  of  a  SIMO  system  as  discussed  above.  The  multiple  transmitted  signals 
Xj  (r,  t ) ,  j  =  1, 2, . . Mt ,  are  convolved  with  the  channel  impulse  responses 

hj  ( T,t ),  j  =  1,2 .  The  receiving  antenna  receives  a  superposition  of  the  multian¬ 
tenna  transmissions  through  the  channels  [12].  Accordingly,  the  received  signal  can  be 
expressed  as 

r (t)  =  xl(t)*hl(T,t)  +  x2(t)*h2(T,t)  +  ...  +  . xM t{t)*hM  (r,t)  .  (2.6) 

Using  the  vector  notation  for  the  received  signal  equation  (2.6),  we  can  write  the  input- 
output  relation  for  a  MISO  system  as 

r(t)  =  h(r,t)*x(t)  (2.7) 

where  x(t)  =  (/)  x2  (/)  •  •  •  xM  (t)J  is  the  transmission  vector,  and  the  channel 

vector  is  denoted  as  h(t)  =  (t)  h2  (t)  •  •  •  hM  (r)J  . 


4.  Multiple-Input  Multiple- Output  (MIMO)  Systems 

After  discussing  the  input-output  relations  for  the  SIMO  and  the  MISO  systems, 
we  now  proceed  to  develop  the  MIMO  system  relations.  A  simple  MIMO  system  is  illus¬ 
trated  in  Figure  4.  The  received  signal  at  the  first  receive  antenna  is  expressed  by 
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(2.8) 


n  (0  =  *1  (0  *  M7>0  +  •  •  •  +  *Mf  (0  *  Km,  {r,t)  . 

This  is  analogous  to  the  MISO  system  input-output  relation.  The  received  signal  at  the 
Mr- th  receive  antenna  is  given  by 

rMr  (0  =  X1  (0  (^"’0  +•••  +  XM,  (0  (^"’0  •  (2.9) 

The  general  input-output  relation  for  a  MIMO  system  in  matrix- vector  notation  is  given 
by 

r(f)  =  H(r,f)*x(f)  (2.10) 

where  x(t)  =  (t)  x2(t)  •••  xM  (r)J  is  the  Mt  xl  transmission  vector, 

r  ~\T 

r(t)=  rx(t)  r2(t)  •••  rM  (t)  is  the  Mr  xl  receive  vector,  and  H(r,t)  is  the 
MrxMt  channel  matrix  given  by 


Figure  4.  Schematic  of  a  MIMO  System  (After  Ref.  2.) 


B.  ALAMOUTI  SCHEME 

The  Alamouti  transmit  diversity  technique  was  proposed  in  1998  [4].  The  tech¬ 
nique  is  generally  referred  to  as  the  Alamouti  scheme  in  the  literature.  Alamouti  scheme 
is  one  of  the  first  space  coding  schemes  developed  for  the  MIMO  systems.  Since  the 
scheme  has  a  simple  transmit  coding  technique  and  a  simple  decoding  implementation,  it 
gained  considerable  interest. 
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The  Alamouti  scheme  was  developed  for  systems  under  flat  fading  conditions, 
i.e.,  over  frequency  independent  channels.  Therefore,  while  discussing  the  scheme,  the 
indices  of  delay  t  and  time  1  of  the  channels  are  dropped.  The  channels  are  denoted  as 
hi ,  i  =  1, 2 ,  representing  single-tap  impulse  responses.  Prior  to  developing  the  Alamouti 

scheme,  we  first  need  to  present  the  input-output  relations  of  the  SIMO  system  using  a 
maximal  ratio  combining  scheme. 

1.  Maximal  Ratio  Combining  (MRC)  Scheme  for  1x2  Systems 

Maximal  ratio  combining  (MRC)  scheme  is  developed  for  the  systems  having 
multiple  receive  antennas,  i.e.,  multiple  channels.  The  scheme  is  based  on  the  assumption 
that  the  receiver  has  perfect  channel  knowledge.  The  transmitted  information  is  estimated 
by  processing  the  channel  state  information  (CSI)  and  the  received  signals  [5].  The  chan¬ 
nel  information  may  be  obtained  by  inserting  known  pilot  symbols.  The  receiver  esti¬ 
mates  the  channel  information  by  interpolating  the  samples  of  the  received  pilot  symbols 

[4]. 


Figure  5  illustrates  a  1x2  system  employing  MRC.  The  impulse  response  of  the 
channel  from  the  transmit  antenna  to  the  receive  antenna  1  is  denoted  by  /y  and  to  the 

receive  antenna  2  is  denoted  by  h2 .  Assuming  that  the  channels  have  flat  fading,  received 
signals  with  additive  noise  are  expressed  by 

rx  =  hlsl  +nl 


r2  =  h2s1  +  n2 


(2.12) 


where  5,  is  the  transmitted  information  symbol  and  n,  and  n2  denote  the  complex  noise 
components.  The  MRC  uses  the  CSI  and  the  received  signals  r{  and  r2  to  compute  the 


estimated  value  of  y .  The  MRC  obtains  an  estimate  of  y  using  the  relation  [4] 

y  =  h* i]  +  h*2r2 .  (2.13) 

The  combining  scheme  compensates  for  the  phase  shift  in  the  channel  by  multiplying  the 
received  signal  with  the  complex-conjugate  of  the  corresponding  channel  [5].  Substitut¬ 
ing  (2.12)  into  (2.13),  we  can  rewrite  the  estimate  ofy  as 

y  =  |2  + 1^2 )  5i  “I-  Kn\  +  Kni  •  (2. 14) 
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The  constructed  complex  value  is  the  maximum  likelihood  estimate  of  the  transmitted 
symbol  [4].  The  estimate  of  the  transmitted  symbol  is  mostly  dependent  upon  the  magni¬ 
tudes  of  the  channels  /r,  and  h2 ,  i.e.,  MRC  is  resistant  to  phase  changes  of  the  channels. 


Figure  5.  Maximal  Ratio  Combining  for  a  1x2  System  (After  Ref.  4.) 

2.  Alamouti  Scheme  for  2x1  Systems 

The  Alamouti  scheme  for  a  system  with  two  transmit  antennas  and  a  single  re¬ 
ceive  antenna  is  shown  in  Figure  6.  Two  consecutive  symbols  5,  and  s2  are  transmitted 
simultaneously  during  the  first  symbol  period  (at  time  =  t).  During  the  next  symbol  pe¬ 
riod  (at  time  =  t  +  T),  -s*2  and  s*  are  transmitted  from  antenna  1  and  2,  respectively, 
where  the  asterisk  *  indicates  complex  conjugation.  The  transmission  sequence  is  shown 


in  Table  1. 


Figure  6.  Alamouti  Scheme  for  a  2x2  System  (After  Ref.  4.) 
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Time 

Transmit  An¬ 
tenna-1 

Transmit  An¬ 
tenna-2 

t 

5i 

S2 

t+T 

* 

-*2 

* 

Table  1.  Transmission  Sequence  of  the  Alamouti  Scheme 


The  received  signal  at  time  t,  which  is  the  superposition  of  the  two  incoming  sig¬ 
nals,  can  be  expressed  by 

t\  =r(t)  =  h1sl  +h2s2  +nl  (2.15) 

where  n]  is  the  additive  complex  noise  component.  Assuming  that  the  channel  is  constant 
across  two  consecutive  symbol  periods,  we  can  express  the  received  signal  at  time  t+T  as 

ri  =  r{t  +  T)  =  h2s*x  —hxs2  +n2  (2.16) 

where  n2  is  the  additive  complex  noise  component. 


Two  consecutive  received  signals  and  the  CSI  are  employed  at  the  Alamouti  de¬ 
coder.  The  estimates  of  the  transmitted  symbols  are  computed  by  [4] 


=  l\rx  +  h2r* 
~s2  =  Krx-hxr*. 


(2.17) 


Substituting  (2.15)  and  (2.16)  into  (2.17),  we  can  rewrite  the  decoder  equations  as 


=  (h  f  +  \K  I' )  +/</?!+  h2n2 

=(M2  +\h2\2)s2+Kni  ~Vv 


(2.18) 


This  outcome  shows  that  the  Alamouti  decoder  compensates  the  phase  changes  of  the 
channel  in  a  manner  similar  to  that  in  the  MRC  scheme. 

3.  Alamouti  Scheme  for  2x2  Systems 

We  can  extend  the  2x1  Alamouti  scheme  to  a  2x2  system  by  adding  one  more  re¬ 
ceiving  antenna.  A  schematic  of  the  Alamouti  scheme  for  a  2x2  system  is  illustrated  in 
Figure  7.  The  transmission  sequence  is  identical  to  the  one  in  the  2x1  system  (see  Ta¬ 
ble  1),  but  the  decoding  scheme  will  be  different. 
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Figure  7.  Alamouti  Scheme  for  a  2x1  System  (After  Ref.  4.) 


Assuming  that  the  channel  is  constant  across  consecutive  symbol  periods,  the  re¬ 
ceived  signals  can  be  expressed  by 

rl=hllsl+hl2s2  +  nl 
r2=h2lsl+h22s2+n2 

(2-lc 

r3  =  ~hns2  +  h1  25,  +n3 

t  —  — II  V*  4-  h  V  *  +  77 

'4  n  2,1J2  '  L  2,2^1 

where  ,  i  =  1, 2, 3, 4 ,  are  the  complex  noise  components.  Table  2  lists  the  notation  for 
the  received  signals  at  the  receive  antennas.  The  channel  impulse  responses 
hs  . ,  i,  j  =  1, 2 ,  represent  the  four  possible  paths  in  the  scheme. 


Time 

Receive 

Antenna- 1 

Receive 

Antenna-2 

t 

ri 

r2 

t+T 

r3 

r4 

Table  2.  Notation  for  the  Received  Signals  in  2x2  Alamouti  Scheme 


The  Alamouti  decoder  constructs  the  estimates  of  v,  and  ,v2  using  the  CSI,  the  re¬ 
ceived  signals  and  the  decoding  equations  given  by  [4] 
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(2.20) 


Clearly,  we  can  see  that  the  estimates  mostly  depend  on  the  magnitudes  of  the  channels 
so  that  the  scheme  is  resistant  to  phase  changes. 

C.  SUMMARY 

This  chapter  discussed  the  input-output  relations  in  MIMO  systems,  the  MRC 
scheme  for  2x1  systems  and  the  Alamouti  scheme  for  2x1  and  2x2  systems.  In  all  cases, 
only  flat  channel  models  are  used  to  introduce  the  decoding  process.  In  the  next  chapter, 
MIMO  channel  models  for  the  frequency  selective  cases  are  introduced. 
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III.  MIMO  CHANNEL  MODELS 


Building  a  channel  model  is  important  for  investigating  the  performance  of  wire¬ 
less  communication  systems.  Channel  parameters  are  extracted  from  measurements  in¬ 
cluding  path  loss,  Ricean  A'- fact  or,  delay  spread  and  Doppler  spread  for  SISO  systems 
and  spatial  correlation  for  MIMO  systems.  The  Stanford  University  Interim  (SUI)  chan¬ 
nels  are  widely  used  channel  models  based  on  such  measurements.  In  the  literature,  most 
simulations,  especially  the  IEEE  802.16a  standard-based  OFDM  system  simulations,  are 
conducted  using  the  SUI  channel  models.  Generation  of  the  SUI  channels  for  SISO  sys¬ 
tems  and  the  spatial  correlation  effects  on  correlated  MIMO  channels  are  discussed  in 
this  chapter. 

A.  STANFORD  UNIVERSITY  INTERIM  (SUI)  CHANNELS  FOR  THE  IEEE 

802.16a  STANDARD 

The  IEEE  802.16  Broadband  Wireless  Access  Working  Group  proposed  channel 
models  for  fixed  broadband  systems  in  the  year  2001.  These  channel  models  are  known 
as  Stanford  University  Interim  (SUI)  channel  models  because  of  the  contribution  of  the 
Stanford  University  in  the  development  process.  The  detailed  channel  generation  is  de¬ 
scribed  in  [11]. 

Three  different  terrain  categories  are  defined  in  the  SUI  channel  models.  The 
maximum  path  loss  category  is  hilly  terrain  with  moderate-to-heavy  tree  densities  (Cate¬ 
gory  A).  The  minimum  path  loss  category  is  typically  flat  terrain  with  light  tree  densities 
(Category  C).  The  intermediate  path  loss  category  is  Category  B.  The  SUI  channel  mod¬ 
els  define  two  models  for  each  terrain  type,  thus  leading  up  to  a  total  of  six  categories. 
Each  channel  model  is  characterized  by  a  time  delay,  a  Ricean  V- fact  or  and  a  Doppler 
spread  in  addition  to  the  path  loss  behavior.  The  broad  parameterization  of  the  SUI  chan¬ 
nels  is  summarized  in  Table  3.  From  this  classification,  we  can  easily  see  that  the  SUI-1 
channel  is  the  least  severe  channel  model  with  a  low  delay  spread  ( 0.4  jus  for  the  2nd  tap 
and  0.9  jus  for  the  3rd  tap),  a  low  Doppler  spread  (0.4  Hz,  0.3  Hz,  0.5  Hz  for  the  1st,  2nd 
and  3rd  taps,  respectively)  and  a  high  V-factor  (8  dB  for  the  1st  tap).  The  most  severe 
channel  is  certainly  the  SUI-6  model  with  a  high  delay  spread  ( 4  /is  for  the  2nd  tap  and 
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10  //s  for  the  3rd  tap),  a  high  Doppler  spread  (2  Hz,  1.5  Hz,  2.5  Hz,  for  the  1st,  2nd  and  3rd 
taps,  respectively)  and  K  =  0 ,  i.e.,  no  direct  path  component  but  only  a  fading  component 
(see  Appendix  A  for  details  of  all  SUI  model  parameters). 


SUI-Model 

Terrain  Type 

Delay  Spread 

Ricean  K  Factor 

Doppler 

SUI-1 

C 

Low 

High 

Low 

SUI-2 

C 

Low 

High 

Low 

SUI-3 

B 

Low 

Low 

Low 

SUI-4 

B 

Moderate 

Low 

High 

SUI-5 

A 

High 

Low 

Low 

SUI-6 

A 

High 

Low 

High 

Table  3.  SUI  Channel  Characterization  (After  Ref.  13.) 


Each  of  the  six  SUI  channels  is  modeled  as  a  3-tap  discrete  system  with  non- 
uniform  delays,  i.e.,  multipath  fading  exists  for  each  SUI  model.  The  SUI  channels  are 
also  modeled  for  different  antenna  types.  Each  channel  has  parameters  defined  for  an 
omni-directional  antenna  and  for  a  30°  directional  antenna.  The  detailed  parameters  for 
the  omni-directional  and  for  the  30°  directional  antenna  case  are  given  in  Appendix  A. 

The  latest  IEEE  802.16  Broadband  Wireless  Access  Working  Group  report  [11] 
contains  a  complimentary  Matlab  code  for  generating  channels  for  the  SUI-3  model.  In 
this  thesis,  to  create  a  simulation  environment,  an  expanded  version  of  Matlab  code  (con¬ 
taining  all  six  SUI  channels)  was  developed  based  on  the  code  in  [11]. 

After  generating  the  SUI  channel  model  in  Matlab,  the  frequency  responses  of  the 
channel  models  were  investigated.  The  frequency  response  plots  of  five  different  SUI 
channels  are  shown  in  Figure  8.  The  frequency  range  was  chosen  to  be  in  terms  of  the 
subcarrier  index  0-255  in  accordance  with  the  fast  Fourier  transform  (FFT)  size  of  the 
IEEE  802.16a  standard.  The  frequency  response  was  normalized  such  that  the  maximum 
magnitude  was  1.  By  observing  the  fluctuations  of  the  magnitude  values,  the  SUI-1 
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channel’s  magnitude  response  changes  smoothly  between  0.9  and  1.0  where  as  the  SUI-5 
channel’s  magnitude  response  varies  rapidly  between  0.5  and  1.0.  As  a  result,  the  higher 
the  SUI  channel  index,  the  more  frequency  selective  is  the  channel. 


Figure  8.  The  Frequency  Response  Plots  of  Different  SUI  Channel  Models 

B.  GENERATING  CORRELATED  MIMO  CHANNELS 

Since  there  are  multiple  antennas  at  both  ends  of  MIMO  systems,  the  presence  of 
spatial  correlation  has  an  effect  on  the  channels.  The  ideal  case  is  that  all  channels  in  a 
MIMO  system  are  independent,  i.e.,  uncorrelated  with  each  other.  In  reality,  however,  the 
physical  configuration  of  the  antennas  (antenna  spacing)  and  operating  environment  pa¬ 
rameters  (angle  of  arrival  and  angle  of  departure  of  the  electromagnetic  waves)  cause 
spatial  correlation  on  the  antenna  array  of  a  MIMO  system.  A  MIMO  system’s  perform¬ 
ance  is  reduced  by  this  spatial  correlation.  An  analysis  of  spatial  correlation  is  discussed 
in  [14]. 

1.  Uncorrelated  MIMO  Channels 

Consider  a  system  with  Mt  transmit  and  Mr  receive  antennas.  Assuming  that  all 
the  channels  are  flat  and  uncorrelated  (i.e.,  each  path  is  represented  by  a  single  complex 
tap  weight),  the  MIMO  channel  can  be  modeled  to  be  zero  mean  and  circularly  symmet- 
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ric  with  unit  variance.  The  uncorrelated  channel  is  an  identically  and  independently  dis¬ 
tributed  channel,  denoted  asHM,  with  dimensions  Mr  xMt .  The  properties  of  Hu,  are 


KH 

{|H,  |}  =  0  if  ^  m  or./  *  n 


where  £{•}  denotes  the  expectation  operator  [2]. 


(3.1) 


2.  Correlated  MIMO  Channels 

In  practice,  MIMO  channels  can  differ  from  Hm  because  of  the  spatial  correla¬ 
tion.  A  more  practical  MIMO  channel  can  be  expressed  as  [2] 

H  =  R‘  H  R12  (3.2) 

where  R,.  is  the  MrxMr  receive  correlation  matrix  and  R  is  the  MtxMt  transmit  cor¬ 


relation  matrix.  The  correlation  matrices  Rr  and  R,  are  positive  semi-definite  Hermitian 
matrices  and  have  the  following  properties: 

R  =R  1/2R  1/2 


R  =R1/2R12. 


(3.3) 


The  correlation  matrices  for  a  2x2  system  can  be  generated  in  the  following 
form  [11]: 

r,.=  1  p;  ,  0.4) 

L  Pr  1 J 

i  p, 

R,  =  7  (3.5) 

La  1 J 

where  pr  and  pt  denote  the  receive  and  transmit  correlation  coefficients.  The  spatial  fad¬ 


ing  coefficient  is  dependent  upon  the  distance  between  the  antennas,  the  wavelength  and 
the  geometry  of  the  physical  environment  and  takes  values  between  0  and  1 .  When  the 
transmit  and  receive  correlation  coefficients  are  both  equal  to  zero,  the  correlation  matri¬ 
ces  R(  and  R;  become  identity  matrices,  and  then  the  MIMO  channel  will  be  equal  to  the 


uncorrelated  channel  Hu, .  The  computation  of  the  correlation  coefficients  is  presented  in 
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detail  in  [14].  In  this  thesis,  the  environmental  parameters  (antenna  spacing,  angle  of  arri¬ 
val  and  angle  of  departure  of  the  electromagnetic  waves)  were  not  taken  into  account  to 
compute  the  correlation  coefficients.  The  correlation  coefficient  parameters  were  chosen 
arbitrarily. 

3.  LOS  Component 

In  the  presence  of  a  line-of-sight  (LOS)  component,  MIMO  channels  have  a  fixed 
and  a  random  component.  The  MIMO  channel  with  a  LOS  component  can  be  represented 
as  [2] 

H 

where  K I {\  +  K)W  is  the  LOS  and  ^j\/(\  +  K)H  w  is  the  random  uncorrelated  fading 

component.  The  elements  of  H  have  unit  power.  The  Ricean  ^-factor  in  (3.6)  is  the  ratio 
of  the  LOS  component  power  to  the  fading  component  power.  When  K  -  0 ,  the  channel 
becomes  pure  Rayleigh  faded  channel,  i.e.,  there  is  no  LOS  component.  When  K  =  ,  all 

we  have  is  the  LOS  component,  and  there  is  no  fading  component.  Then  the  channel  be¬ 
comes  flat  [2] . 

In  the  presence  of  both  spatial  correlation  and  LOS  component,  to  achieve  a  com¬ 
plete  MIMO  channel,  only  the  fading  component  is  computed  using  (3.2),  and  the  fixed 
LOS  component  is  added  to  the  fading  component  using  (3.6). 

The  generation  of  the  correlated  SIMO  and  MISO  channels  are  similar  to  the 
MIMO  channels.  Consider  a  system  with  a  single  transmit  and  two  receive  antennas.  The 
uncorrelated  SIMO  channel  Hu,  is  a  2x1  row  vector.  Since  we  have  a  single  transmit  an¬ 
tenna,  there  is  no  transmit  correlation,  i.e.,  pi  =  0 .  The  transmit  covariance  matrix  R,  is 
then  an  identity  matrix.  From  (3.2),  the  correlated  SIMO  channel  can  be  expressed  as 

H  =  R/2Hh  (3.7) 

where  R1/2  is  the  square-root  of  the  receive  correlation  matrix  R  . 

For  the  MISO  case,  consider  a  system  with  two  transmit  antennas  and  a  single  re¬ 
ceive  antenna.  The  uncorrelated  MISO  channel  H  will  be  a  1x2  column  vector.  Since 

w 
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there  are  multiple  antennas  only  at  the  transmitter,  the  receive  correlation  matrix  Rr  is  an 
identity  matrix.  The  correlated  MISO  channel  is  expressed  by 

H  =  HmR|2  (3-8) 

where  R'/2  is  the  square-root  of  the  transmit  correlation  matrix  R  . 


In  the  presence  of  a  LOS  component,  the  fixed  components  of  the  MISO  and  the 
SIMO  channel  are  added  to  the  correlated  fading  components  of  the  MISO  and  the  SIMO 
channels,  respectively.  The  complete  MISO  or  SIMO  channel  is  achieved  by  the  relation 
given  in  (3.6). 

C.  GENERATING  CORRELATED  FREQUENCY-SELECTIVE  MIMO 

CHANNELS 

The  generation  of  the  SUI  channels  for  SISO  systems  and  correlated  flat  MIMO 
channels  have  been  discussed  so  far.  The  uncorrelated  channel  Hh  in  (3.2)  has  single-tap 

complex  components  for  the  flat  channel  case.  The  frequency  selective  channels  have 
multi-tap  components.  To  study  the  correlated  frequency  selective  MIMO  channels,  con¬ 
sider  a  system  with  two  transmit  and  two  receive  antennas  operating  over  the  MIMO 
channels  generated  by  the  SUI  models.  Each  tap  is  characterized  by  a  delay  spread,  a 
power  level,  a  Doppler  spread  and  a  Ricean  V-factor.  Therefore,  a  SUI  channel  has  a 
fixed  and  a  fading  component  for  each  tap.  The  fading  component  includes  the  delay 
spread,  the  Doppler  spread  and  the  power  level.  The  fixed  component  is  generated  with 
the  specified  Ricean  V- factor  for  the  given  SUI  model  [1 1],  In  a  MIMO  channel  for  a  2x2 
system  based  on  SUI  channel  modeling,  the  fading  and  fixed  components  can  be  ex¬ 
pressed  as 


"MO  M2)  M3)  MO  M2)  M3 )" 
.MO  M2)  *2,(3)  A,,  (0  *22(2)  *22(3). 


H  = 


MO  M2)  M3) 

*2,(0  *2,  (2)  *2,(3) 


*, 2(1)  M2)  *12(3) 

*22(0  *22(2)  *22(3).’ 


(3.9) 

(3.10) 


respectively,  where  h:j  (k  )  is  the  channel  response  between  the  /-th(z  =  1,2)  transmit  an¬ 
tenna  and  the  j-th  ( j  =  1, 2)  receive  antenna  for  the  A-th  ( A  =  1, 2, 3)  tap.  Assuming  that  all 
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tap  correlations  are  equal,  the  correlated  fading  component  of  the  MIMO  channel  is  ex¬ 
pressed  by 

hs(*)=r:,2h„,(or:'2  (3.ii) 

where  Hu  (A)  is  the  2x2  uncorrelated  channel  matrix  for  the  k- th  tap. 

To  achieve  the  complete  MIMO  channel,  the  LOS  component  will  be  added  to  the 
correlated  fading  component.  Then  the  MIMO  channel  is  given  by 

=  (3'12) 

where  H;/  (A)  is  the  fixed  component  between  the  /-th  (/  =  1,2)  transmit  antenna  and  the 
y-th  (7  =  1, 2)  receive  antenna  for  the  A-th  ( A  =  1, 2, 3)  tap.  Generation  of  a  correlated 

MIMO  channel  is  illustrated  in  Figure  9. 

D.  SUMMARY 

In  this  chapter,  the  generation  of  SUI  channels  for  SISO  systems,  uncorrelated 
MIMO  channels  and  spatial  correlation  effects  on  flat  and  frequency  selective  MIMO 
channels  were  discussed.  It  was  demonstrated  that  the  channel  becomes  more  frequency 
selective  as  the  SUI  index  increases. 

In  the  next  chapter,  we  will  develop  the  input-output  relations  for  the  Alamouti 
scheme  based  STBC  systems  over  frequency  selective  multi-tap  channels. 
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Power  level 


Figure  9.  Conceptual  Diagram  for  Generating  a  Correlated  MIMO  Channel  Based  on 

SUI  Channel  Models 
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IV.  ALAMOUTI-BASED  SCHEMES  OVER  FREQUENCY- 

SELECTIVE  CHANNELS 

The  transmitted  signals  experience  different  types  of  fading  depending  upon  the 
relationship  between  the  signal  parameters  (bandwidth,  symbol  period,  etc.)  and  the 
channel  parameters  (Doppler  spread,  delay  spread,  etc.).  When  the  bandwidth  of  the 
channel  is  greater  than  the  bandwidth  of  the  transmitted  signal,  the  received  signal  ex¬ 
periences  a  flat  fading  channel.  A  flat  channel  is  approximated  as  a  single-tap  weight 
with  zero  delay.  When  the  bandwidth  of  the  channel  is  greater  than  the  bandwidth  of  the 
transmitted  signal,  the  received  signal  experiences  a  frequency  selective  channel  or  a 
multipath  fading  channel.  Over  frequency  selective  channels,  multiple  versions  of  the 
transmitted  signal  having  faded  amplitudes  arrive  at  the  receiver  at  slightly  different  time 
instants.  Frequency  selective  channels  are  modeled  by  multiple  tap  discrete  systems  with 
different  delay  parameters  [15]. 

Alamouti’s  transmit  diversity  scheme,  which  was  discussed  in  Chapter  II,  was 
originally  proposed  for  flat  channels.  The  space  coding  technique  is  performed  over  a 
symbol  pair,  i.e.,  two  consecutive  symbols  (see  Table  1).  The  decoding  equations  (2.17) 
and  (2.20)  are  developed  for  flat  channels.  The  frequency-flat  channels  in  these  equations 
are  represented  as  one-tap  complex  weights.  The  extension  of  the  Alamouti  scheme  to 
frequency- selective  channels  was  investigated  in  [8]  and  [9].  The  extended  scheme  is 
called  space-time  block  coding  (STBC),  since  space  coding  is  performed  over  blocks  of 
symbols.  The  scheme  presented  in  [8]  and  [9]  is  for  the  systems  with  two  transmit  anten¬ 
nas  and  a  single  receive  antenna  (2x1).  In  this  thesis,  the  scheme  is  further  extended  to 
1x2  and  2x2  systems  using  the  original  Alamouti  scheme. 

OFDM  has  emerged  as  an  attractive  alternative  scheme  to  conventional  modula¬ 
tion  schemes  due  to  its  effectiveness  in  reducing  the  effects  of  frequency  selective  chan¬ 
nels  [7].  The  combination  of  MIMO  systems  and  OFDM  modulation  is  one  of  the  widely 
discussed  areas  in  the  MIMO  research.  Different  space  coding  techniques  in  MIMO- 
OFDM  systems  have  been  investigated  by  many  authors  [6,  7,  16].  In  this  work,  we  have 
selected  the  STBC  scheme  as  the  space  coding  scheme  for  MIMO-OFDM  systems  to 
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study  its  simulation  performance  and  compare  its  performance  to  that  of  STBC-SC  sys¬ 
tems.  The  STBC  scheme  is  chosen  as  an  optional  technique  in  the  OFDM  physical  layer 
specifications  of  the  IEEE  802.16a  standard. 

In  this  chapter,  we  discuss  the  STBC  coding  technique  for  the  MIMO  systems 

with  both  single-carrier  modulation  and  OFDM.  Input-output  relations  and  the  decoding 

equations  are  presented  for  SISO,  SIMO,  MISO  and  MIMO  systems. 

A.  SPACE-TIME  BLOCK  CODING-SINGLE  CARRIER  (STBC-SC)  SYS¬ 
TEMS  OYER  FREQUENCY-SELECTIVE  CHANNELS 

Before  discussing  the  STBC  schemes,  in  order  to  get  familiar  with  the  notation, 

the  discrete-time  SISO  system  relations  over  frequency  selective  channels  are  introduced 

first. 

1.  SISO  Systems  over  Frequency-Selective  Channels 

Consider  that  a  block  of  symbols  is  transmitted  from  a  single  antenna  by  assum¬ 
ing  that  the  channel  characteristics  do  not  change  during  the  period  of  transmission.  The 
transmitted  block  of  symbols  is  expressed  by 

S  =  [s(l)  ^(2)  ...  5(A)]  (4.1) 

where  N  is  the  block  size.  The  transmitted  symbol  block  will  experience  a  frequency  se¬ 
lective  channel.  The  channel  can  be  represented  as  a  discrete-time  filter  given  by 

h  (q-1 )  =  h0  +  \q-1  + . . .  hPq~P  (4.2) 

where  q  1  is  the  unit  delay  operator  and  hn  ,n  =  l,2,...P  are  the  filter  coefficients.  The 
received  signal  can  be  expressed  as 

r(t)  =  h (<7  1 )  v (?)  =  h0s(t)  +  h{s(t -1)  + ...hps(t - P) .  (4.3) 

The  noise  component  of  the  received  signal  is  omitted  here  for  convenience.  The  maxi¬ 
mum  likelihood  estimate  of  the  transmitted  symbols  block  is  computed  by 

s  (t)  =  h*  ( q)r(t )  =  h*r(t)  +  h^r(t  +  l)  +  ...  +  hpr(t  +  P)  (4.4) 

where  h*  (q)  represents  a  non-causal  realization  of  the  filter  given  in  (4.2)  [8]. 

2.  SIMO  Systems  over  Frequency- Selective  Channels 

The  schematic  of  a  SIMO  system  over  frequency  selective  channels  is  illustrated 
in  Figure  10  for  a  1x2  realization.  The  transmitted  block  of  symbols  is  identical  to  the 
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SISO  case  (see  Eq.  (4.1)).  Since  there  are  two  receive  antennas,  we  have  two  channel  fil¬ 
ters  given  by 


\(q~l)  =  \  +\q~'  +---+\Pq~p 
h~i  ( q  ) —  ^20  ^2  q  +  •  •  • + hlp  q 


(4.5) 


The  transmitted  signal  is  passed  through  these  two  channel  filters,  and  the  received  sig¬ 
nals  are  expressed  by 


rl(t)  =  hl(q~1)s(t) 

(4.6) 

r2(t)  =  h2(ql)s(t). 


Channel 

State 

Information 


Figure  10.  Schematic  of  a  1x2  SIMO-SC  System  Utilizing  MRC  over  Frequency- 

Selective  Channels 


The  estimate  of  the  transmitted  signal  block  is  computed  on  the  lines  of  the 
maximal  ratio  combining  scheme.  The  multiplication  operation  in  the  decoding  equation 
(2.13)  is  turned  into  a  filter  operation  for  the  frequency  selective  channels.  Accordingly, 
the  estimate  of  the  transmitted  symbol  block  is  given  by  [8] 

*{t)  =  /<  (q)rt  ( t)  +  h*  ( q)r2  ( t ) .  (4.7) 

Notice  that  the  conjugates  of  the  filters  l\  (q)  and  ti2  (q)  in  (4.7)  are  non-casual  filters. 

By  substituting  (4.5)  into  (4.7),  the  estimate  of  the  transmitted  signal  block  in  (4.7)  can  be 
rewritten  as 

s  (0 = Kr i  (?) + K  6  (? + !) + •  •  • + K/x  {t + p) + K/2  (0  +  (4  g) 

h*r2(t  +  l)  +  ...  +  h*2pr2(t  +  p). 
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3.  MISO  Systems  over  Frequency- Selective  Channels 

The  Alamouti  scheme  for  a  2x1  system  over  flat  channels  was  discussed  in  Chap¬ 
ter  II.  The  flat-channel  implementation  was  based  on  applying  the  Alamouti  encoding 
scheme  over  a  symbol  pair,  i.e.,  two  consecutive  symbols.  The  scheme  can  be  extended 
to  the  frequency- selective  channel  case  by  utilizing  the  STBC  technique.  A  schematic  of 
a  2x1  MISO  system  utilizing  STBC  is  illustrated  in  Figure  11.  The  transmitted  blocks  of 
symbols  are  given  by 

S,=p,(l)  j,(2)  ...  S,(Ar)] 

S,  =[.*,(!)  (2)  ...  (AT)] 

where  S,  and  S2  are  the  simultaneously  transmitted  blocks  of  symbols  from  antenna  1 

and  antenna  2,  respectively.  The  Alamouti  space-time  encoding  scheme  is  performed 
over  these  two  symbol  blocks.  The  first  antenna  transmits  a  complex-conjugated  and 
symbol-inverted  version  of  S2,  and  the  second  antenna  transmits  a  complex-conjugated 

version  of  S,  in  the  second  transmission  burst.  The  transmitted  blocks  of  symbols  in  the 
second  burst  are  given  by 

s,  =-S;  =[-■!,*  (1)  -.<(2)  ...  -,<(A0] 
s4=s;=[>;(i)  4(2)  ...  «w)] 

where  S3  and  S4are  transmitted  from  antenna  1  and  2,  respectively.  The  transmission 
symbol  block  sequence  for  the  STBC  for  a  2x1  system  is  shown  in  Table  4. 


Figure  11.  Schematic  of  a  2x1  STBC-SC  System  over  Frequency-Selective  Channels 
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Time 

Transmit  An¬ 
tenna-1 

Transmit  An¬ 
tenna-2 

1st  transmission  burst 

s2 

2nd  transmission  burst 

S3  =  -s; 

s4=s; 

Table  4.  Transmission  Sequence  for  STBC 


The  notation  of  the  channels  for  a  2x1  system  is  identical  to  the  ones  for  the  1x2 
system  (see  Eq.  (4.5)).  Assuming  that  the  channel  characteristics  remain  constant  during 
consecutive  transmissions  of  symbol  blocks,  the  received  signals  are  given  by  [8] 

n  (0  =  \  (q~l )  Si  (0 + K  (cf' )  s2  (t) 

(4.11 

r2  (0  =  K  (Cf'  )  S3  {t)  +  h2  ( cf 1  )  54  (t) 

where  r{  (/)  and  r2  (/)  denote  the  received  signals  during  the  1st  and  2nd  transmission 


bursts,  respectively.  The  2x1  STBC  decoder  estimates  the  transmitted  blocks  of  symbols 
using  the  decoding  equations  given  by  [4,  8] 


Sl  (0  =  K  U)ri  (0  +  /E  (Cf')r2  (0 

s~2 (0  =  K  {q)ri  (0 (cf')ri  (0 


where  (q)  and  K(q)  represent  the  non-causal  realizations  of  the  channel  filters. 

4.  MIMO  Systems  over  Frequency- Selective  Channels 

After  discussing  the  SISO,  SIMO  and  MISO  systems  over  frequency- selective 
channels,  we  now  proceed  to  develop  the  MIMO  system  relations.  A  schematic  of  a  2x2 
MIMO  system  utilizing  STBC  is  illustrated  in  Figure  12.  The  spatial  encoding  scheme  is 
identical  to  the  MISO  system  utilizing  STBC.  The  transmitted  blocks  of  symbols  were 
introduced  in  Equations  (4.9)  and  (4.10).  The  four  MIMO  channels  for  the  2x2  system 
are  expressed  as 
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(4.13) 


\,i  (#  ) —  \i0  ^i.i, ^  +...  +  /^lp^ 

\,2  (*7  )  —  ^l,20  "*"  ^1.2,^  +  . .  . +  /Z[  2(J<7 
/l2  2  («")  ^2,20  +  K,2tCl  '  +  -'-  +  ^2,2p<? 
^2,2  ^2,20  +  K,2xC1  '  +  -..  +  ^2,2p<? 


Figure  12.  Schematic  of  a  2x2  STBC-SC  System  over  Frequency-Selective  Channels 


Assuming  that  the  channel  characteristics  remain  constant  over  two  transmission  bursts, 
the  received  signals  are  given  by 

ri  (0  =  hu  W' )  (0  +  Ki  [<fl )  *2  (0 

r2  (0  =  K\  (<T‘ )  Si  (0  +  K2  W )  s2  (0 

(4.F 

r3  {<)  =  Ki  Wx )  53  (0 + Ki  (cf' )  -c  (0 

r4  (0  =  Kl  (Cl  l  )  S3  (0  +  K,2  (<7^  )  ^4  (O' 

The  sequence  and  the  location  of  the  received  symbol  blocks  in  (4.14)  are  given  in  Ta¬ 
ble  5.  The  STBC  decoder  estimates  the  transmitted  blocks  of  symbols  by  using  the  de¬ 
coding  equations  [4,  8] 

Si  (0  =  h  *u(q)  #;  (t)  +  h  l2{q~l )  r*  (t)  +  h21  ( q)r2  ( t )  +  h  2  2(q~l )  r*  ( t ) 

(4.1: 

s2  (0  =  KM)ri  (0  ~h  u(q~x  )c  (0  +  h  \M)r2  (0 -h2,i((fl  )K  (0- 


where  h*  j(q),  i  =  1, 2  ;  j  =  1, 2 ,  represent  the  non-causal  realizations  of  the  channel  fil¬ 
ters.  The  equations  in  (4.15)  are  similar  to  the  flat  channel  Alamouti  decoding  equations 
previously  introduced  in  (2.20).  The  multiplication  operation  in  the  flat  channel  case  is 
now  replaced  by  the  filter  operation. 
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Receive  An¬ 
tenna-1 

Receive  An¬ 
tenna-2 

1st  transmission  burst 

n(t) 

2nd  transmission  burst 

rA{) 

Table  5.  The  Sequence  and  the  Location  of  the  Received  Signals  for  a  2x2  STBC 

System 

B.  SPACE-TIME  BLOCK  CODING-ORTHOGONAL  FREQUENCY  DIVI¬ 
SION  MULTIPLEXING  (STBC-OFDM)  SYSTEMS  OVER  FREQUENCY- 
SELECTIVE  CHANNELS 

In  the  previous  section,  we  discussed  the  single-carrier  MIMO  systems  utilizing 
STBC  over  frequency- selective  channels.  The  combination  of  STBC  and  OFDM  is  dis¬ 
cussed  in  this  and  the  following  sections.  OFDM  has  emerged  as  an  attractive  and  alter¬ 
native  scheme  to  the  conventional  modulation  schemes  due  to  its  effectiveness  in  reduc¬ 
ing  the  effects  of  frequency- selective  channels  [7].  In  OFDM,  the  entire  signal  bandwidth 
is  divided  into  a  number  of  orthogonal  subcarrier  bands,  and  then  the  signal  is  transmitted 
in  these  narrowbands  over  a  number  of  subcarriers.  The  OFDM  and  STBC  parameters 
were  selected  from  the  IEEE  802.16a  standard  [10].  The  IEEE  802.16a  standard  includes 
the  2x1  STBC  scheme  as  an  optional  specification  for  systems  employing  OFDM. 

1.  SISO-OFDM  Systems  over  Frequency- Selective  Channels 
A  block  diagram  of  a  SISO-OFDM  system  is  shown  in  Figure  13.  A  detailed  ex¬ 
planation  of  each  block  is  given  in  the  following  discussion.  The  physical  layer  part  of 
the  IEEE  802.16a  standard  is  used  in  describing  the  SISO-OFDM  scheme  here. 
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message)  (estimation 

process) 

Figure  13.  Block  Diagram  of  a  SISO-OFDM  Scheme 

a.  Channel  Encoding 

The  information  bit  sequence  is  first  encoded  by  the  channel  encoder.  Er¬ 
ror  control  coding  is  used  to  rearrange  the  transmitted  information  to  increase  its  resis¬ 
tance  to  noise.  The  IEEE  802.16a  standard  uses  concatenated  forward  error  correction 
(FEC),  which  is  based  on  the  serial  concatenation  of  a  Reed-Solomon  outer  code  and  a 
rate  compatible  Trellis  coded  modulation  (TCM)  inner  code.  Between  the  outer  and  the 
inner  code,  employing  an  interleaver  is  optional.  Since  our  main  focus  in  this  thesis  was 
to  study  the  MIMO  and  MIMO-OFDM  systems,  the  simulation  only  used  the  convolu¬ 
tional  encoder  to  keep  it  simple.  The  encoder  had  a  rate  of  V2,  constraint  length  7  and 
generator  polynomials  in  octal  form  171  and  133  [10]. 

b.  Symbol  Mapping 

The  channel  encoded  bits  are  mapped  to  I  and  Q  symbol  coordinates  using 
a  gray  code  symbol  map  depending  on  the  modulation  scheme.  Figure  14  shows  some  of 

the  constellations  introduced  in  IEEE  802.16a  standard. 
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Figure  14.  Gray  coded  (a)  BPSK,  (b)  QPSK  and  (c)  16-QAM  constellations  (From  Ref. 

10.) 


c.  Inverse  Fast  Fourier  Transform  (IFFT)  Operation 
Mapped  symbols  are  next  input  to  the  IFFT  block.  The  FFT  size  is  equal 
to  the  number  of  subcarriers.  In  the  IEEE  802.16a  standard,  two  types  of  OFDM  schemes 
are  included,  OFDM  with  256  subcarriers  and  OFDMA  (for  higher  data  rates  and  multi¬ 
ple  access)  with  2048  subcarriers.  In  this  thesis,  the  OFDM  with  256  subcarriers  was  cho¬ 
sen  for  simulation  purposes.  Three  types  of  subcarriers  are  defined  in  the  OFDM  tech¬ 
nique  of  IEEE  802.16a  standard:  data  carriers  for  information  transmission,  pilot  carriers 
designed  to  extract  the  channel  information  at  the  receiver,  and  guard  carriers  (also  called 
null  carriers)  placed  on  the  either  edge  of  the  spectrum  to  avoid  interference  from  adja¬ 
cent  bands.  The  assignment  of  subcarriers  in  the  IEEE  802.16a  standard  is  given  in  Table 
6.  These  carrier  allocations  will  be  used  in  organizing  the  information  and  pilot  symbols 
before  the  IFFT  operation.  Figure  15  illustrates  the  subcarrier  organization  prior  to  IFFT 
operation  according  to  the  carrier  allocation  given  in  Table  6. 
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Size  of  FFT 

256 

#  of  information  subcarriers 

192 

#  of  pilot  subcarriers 

8 

#  of  null  subcarriers  (including  the  DC 

subcarrier) 

56 

#  of  lower  frequency  guard  subcarriers 

28 

#  of  higher  frequency  guard  subcarriers 

27 

Frequency  indices  of  null  subcarriers  (in¬ 
cluding  the  DC  subcarrier) 

-128,-127,. ..,-101,0 

+101, +102, ...,+127 

Frequency  indices  of  pilot  subcarriers 

-84,  -60,  -36,  -1 2, 1 2, 36, 60, 84 

Table  6.  Allocation  of  OFDM  Subcarriers  in  the  IEEE  802.16a  Standard 
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Figure  15.  Allocation  of  Subcarriers  before  the  IFFT  Block  (see  Table  6) 


The  output  of  the  IFFT  block  in  discrete  time  can  be  expressed  by 

N  fpt  — 1 


:[n\- 


1 


FFT 


I 

k= 0 


jl!Cfkn 


n  =  0,l...NFFr  -1 


(4.16) 
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where  sk  is  the  A-th  symbol  entering  to  the  IFFT  block  and  fk  is  the  A-th  subcarrier  fre¬ 
quency.  The  index  of  subcarrier  frequency  in  (4.16)  is  chosen  to  be  0  to  NFFT  -1  instead 
of  -128  to  +128 ,  for  convenience.  The  length  of  the  IFFT  output  is  called  the  useful 
time  Tb . 

d.  Cyclic  Prefix  Addition 

Multipath  delay  spread  causes  intersymbol  interference  (ISI),  and  the  ISI 
brings  about  performance  degradation.  To  deal  with  this  problem,  a  guard  time  is  intro¬ 
duced  for  each  OFDM  symbol.  The  guard  time  is  chosen  to  be  larger  than  the  expected 
delay  spread  so  that  multipath  components  from  one  symbol  would  not  interfere  with 
those  from  the  next  symbol.  In  the  IEEE  802.16a  standard,  optional  guard  lengths  Tg  are 

specified  with  respect  to  the  useful  time  7^  .  The  optional  guard  ratios  Tg  jTh  are 

1/4,1/8,1/16  and  1/32  .  The  desired  length  of  the  guard  interval  is  derived  from  this  ratio. 
The  addition  of  cyclic  prefix  (CP)  is  shown  in  Figure  16.  Consider  that  the  guard  ratio  is 
chosen  to  be  1/16  .  The  output  of  the  IFFT  block  has  a  length  of  256;  therefore,  the  guard 
length  will  be  16.  The  CP  is  formed  by  taking  the  last  16  samples  of  the  IFFT  output  and 
concatenating  them  to  the  beginning  of  the  symbol  sequence  in  (4.16).  The  resulting  new 
sequence  is 

Xcp=[x[ 240]  jc[241]  ...  *[254]  *[255]  *[0]  *[l]  ...  *[255]].  (4.17) 
The  total  interval  after  adding  the  CP  is  referred  to  as  the  symbol  time  7  . 


Figure  16.  OFDM  Symbol  Time  Structure  (From  Ref.  10.) 
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e.  Digital-to -Analog  Conversion  and  RF  Modulation 

After  CP  addition,  the  discrete  time  domain  samples  are  converted  to  ana¬ 
log  signals  using  a  digital-to-analog  (D/A)  converter.  The  converted  signal  is  a  continu¬ 
ous  time  baseband  OFDM  symbol.  Then  the  OFDM  symbols  are  upconverted  to  a  radio 
frequency  carrier,  amplified  and  transmitted  through  the  antenna.  The  IEEE  802.16a  stan¬ 
dard  is  specified  for  the  2-11  GHz  band.  In  this  thesis,  the  simulations  were  conducted  at 
the  baseband  level,  i.e.,  no  RF  modulation  was  implemented. 

/.  RF  Demodulation  and  Digital  to  Analog  Conversion 
The  transmitted  OFDM  symbols  experience  the  channel  and  reach  the  re¬ 
ceiver.  Assuming  that  the  channel  is  modeled  as  a  multitap  filter,  the  transmitted  signal  is 
passed  through  a  discrete  time  filter.  The  coefficients  of  the  filter  are  the  complex  tap 
values  of  the  channel.  The  received  signal  now  can  be  expressed  by 

r[t)  =  h(q~1)x(t)  (4.18) 

where  7?  (<7  1 )  is  the  discrete-time  filter  representation  of  the  channel  and  x(t)  is  the 

baseband  OFDM  symbol  transmitted  from  the  antenna. 

g.  Cyclic  Prefix  Removal  and  FFT  Operation 

The  cyclic  prefix  added  at  the  transmitter  is  removed  to  process  the  needed 
information.  After  removal  of  the  CP,  we  have  the  256  samples  carrying  the  transmitted 
symbol  and  channel  information.  These  256  samples  are  processed  through  the  FFT 
block,  and  the  information  symbols  and  the  channel  information  are  obtained.  The  esti¬ 
mate  of  the  symbol  sent  on  the  k- th  subcarrier  is  given  by 

s(k)  =  H*(k)r(k)  (4.19) 

where  H  (k)  is  the  frequency  response  of  the  channel  at  the  k-th  subcarrier  frequency 
and  r(k)  is  k-th  element  of  the  sequence  obtained  by  the  FFT  operation.  In  all  simula¬ 
tions,  we  assumed  that  there  was  full  channel  knowledge  at  the  receiver,  i.e.,  the  channel 
estimation  was  not  performed.  The  frequency  response  of  the  channel  H  (k)  was  ob¬ 
tained  by  taking  the  discrete  Fourier  transform  of  the  channel  impulse  response  with  an 
FFT  of  size  256. 
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In  practice,  channel  state  information  is  carried  to  the  receiver  in  the  pilot 
carriers.  Since  the  channel  estimation  was  not  included  in  our  simulations,  the  pilot  sym¬ 
bols  were  set  to  zero.  Hence  pilot  insertion  was  not  carried  out  in  the  simulation  steps  but 
left  open  for  future  work. 

h.  Demapping  and  Decoding 

The  estimates  of  the  transmitted  symbols  are  then  demapped  to  bits  with 
reference  to  the  constellation  and  sent  to  the  Viterbi  decoder.  The  Viterbi  decoder  is  used 
to  decode  bit  streams  encoded  by  the  convolutional  encoder.  A  Viterbi  decoder  searches 
the  trellis  to  determine  the  path  that  can  create  the  bit  stream  closest  to  the  received  one. 
The  measure  of  closeness  can  be  measured  using  hard  or  soft  decision  decoding  tech¬ 
niques.  The  hard  decision  decoding  technique  tries  to  find  the  bit  stream  path  that  has  the 
minimum  Hamming  distance  from  the  received  sequence.  In  soft  decision  decoding,  the 
bit  stream  path  that  has  the  minimum  Euclidian  distance  from  the  received  sequence  is 
computed  [17].  In  the  simulation  process,  hard  decision  decoding  was  employed  at  the 
Viterbi  decoder.  The  details  of  the  Viterbi  decoder  can  be  found  in  [18]. 

2.  SIMO-OFDM  Systems  over  Frequency-Selective  Channels 
Since  we  have  given  a  detailed  explanation  for  SISO-OFDM  systems  in  the  pre¬ 
ceding  discussion,  we  will  skip  the  following  steps  in  discussing  the  MIMO-OFDM  sys¬ 
tems:  channel  encoding,  symbol  mapping,  pilot  insertion,  digital-to-analog  conversion, 
and  RF  modulation  at  the  transmitter;  RF  demodulation,  analog-to-digital  conversion, 
cyclic  prefix,  extracting  the  information  symbols  from  the  FFT  block,  and  Viterbi  decod¬ 
ing  at  the  receiver.  The  resulting  schematic  of  a  system  utilizing  OFDM  with  single 
transmit  antenna  and  two  receive  antennas  is  illustrated  in  Figure  17. 
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Figure  17.  Schematic  of  a  1x2  SIMO-OFDM  System  Utilizing  MRC 


The  transmitted  OFDM  symbol  will  experience  two  different  channels  1\  and  h2 . 
The  received  signals  are  given  by 

r\  {t)  =  hl(q~1)x(t) 

(4.20) 

r2(t)  =  h2(q~1)x(t) 

where  f\  [q  1  j  and  h2  (</  1  j  are  the  discrete-time  filter  representations  of  the  two  SIMO 

channels  and  x(t)is  the  transmitted  OFDM  signal.  At  the  receiver,  the  CP  is  removed 

and  the  FFT  of  the  256  time-domain  samples  is  taken.  By  using  the  MRC  scheme,  the 
estimate  of  the  transmitted  symbol  at  the  k- th  subcarrier  is  computed  by 

s(k)  =  H;(k)rl(k)  +  H*2(k)r2(k)  (4.21) 

where  H*  (k  )  and  H2  (k  )  are  the  complex  conjugate  versions  of  the  frequency  responses 
of  the  channels  at  the  A-th  subcarrier,  and  r1  ( A )  and  r2  (A)  are  the  A-th  output  of  the  FFT 
blocks. 

3.  MISO-OFDM  Systems  over  Frequency-Selective  Channels 

A  schematic  of  a  MISO-OFDM  system  with  two  transmit  antennas  and  a  single 
receiving  antenna  utilizing  STBC  is  shown  in  Figure  18.  In  SISO-OFDM  and  SIMO- 
OFDM  systems,  only  one  OFDM  symbol  (carrying  192  information  symbols  as  specified 
in  the  IEEE  802.16a  standard)  is  transmitted  through  a  given  transmit  antenna  at  a  time. 
Since  there  are  two  transmit  antennas  for  the  MISO-OFDM,  two  different  OFDM  sym- 
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bols  (making  it  a  total  of  384  information  symbols)  are  transmitted  from  the  two  transmit 
antennas  simultaneously.  The  STBC  encoder  performs  encoding  over  two  consecutive 
symbol  blocks.  The  four  encoded  symbol  blocks  for  transmission  are  constructed  as  fol¬ 
lows: 

5i=[5[!]  s[ 2]  •••  4192]] 

5,  =[5  [193]  j  [194]  ...  5  [384]] 

r  n  (4.22' 

s3=[-/[193]  -5*  [194]  ...  -s*  [384]] 

54=[/[l]  s*  [2]  ...  .v  1 192]]. 

This  block  construction  is  analogous  to  the  STBC-SC  scheme  in  Section  A;  therefore,  in 
the  literature  this  scheme  is  also  called  STBC-OFDM.  The  allocation  of  subcarriers  is 
carried  out  as  specified  in  the  IEEE  802.16a  standard  for  OFDM  (see  Figure  15). 
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Figure  18.  Schematic  of  a  2x1  MISO-OFDM  System  Utilizing  STBC 


Each  of  the  OFDM  symbols  xt  ( t ),  i  =  1,2, 3, 4,  contains  a  block  of  192-symbol, 


sn  i  =  1,2, 3,4  ,  presented  in  (4.22),  respectively.  Assuming  that  the  channel  characteris- 
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(4.23) 


(4.24) 


tics  remain  constant  over  two  consecutive  OFDM  symbols,  the  received  signals  are  ex¬ 
pressed  by 

n  (0  =  \  [t)  +  h2  (q  ')x2  (t) 

r2(t)  =  hl(q~l)x3(t)  +  h2(q~l)x4(t) 

where  l\  (q1 )  and  h2  (q' )  are  the  discrete-time  filter  representations  of  the  two  MISO 

channels.  The  STBC  decoder  at  the  receiver  computes  the  estimates  of  the  transmitted 
signal  at  the  A-th  subcarrier  frequency  using  the  decoding  equations  given  by  [4,  7] 

sl(k)  =  H*  (k)rl(k)  +  H2(k)r*  ( k ) 
s2(k)  =  H*(k)rl(k)~Hl(k)r*  ( k ) 

where  H*  (A)  and  H2  (A)  are  the  complex  conjugates  of  the  frequency  responses  of  the 
channels  at  the  k- th  subcarrier  frequency,  and  /)  ( A )  and  r2  (A)  are  the  A-th  output  of  the 
FFT  blocks. 

4.  MIMO-OFDM  Systems  over  Frequency- Selective  Channels 

A  schematic  of  a  MIMO-OFDM  system  utilizing  STBC  with  two  transmit  and 
two  receive  antennas  is  shown  in  Figure  19.  The  transmitter  side  of  the  system  is  identical 
to  the  MISO-OFDM  employing  STBC.  Assuming  that  the  channel  characteristics  remain 
constant  over  the  period  of  two  transmission  bursts,  the  received  signals  are  given  by 

n  (0  =  Ki  W' )  t  (0 + hi,i  (K )  *2  (0 

r2  (0  =  Kl  W')X  1  it)  +  h2,2  W')X2  (0 
r3  (0  =  Kl  iCf'  )  X3  (0  +  h2A  (Cf'  )  X4  (0 
r4  (0  =  Kl  )  X3  (0  +  h2,2  (K  )  X4  (0 

where  k  t(  q  1 ) ,  f  =  1,2  j  =  1, 2 ,  are  the  discrete-time  filter  representations  of  the  channel 

from  the  j- th  transmit  antenna  to  the  i-th  the  receive  antenna.  The  STBC  decoder  at  the 
receiver  antenna  will  estimate  the  transmitted  symbols  at  the  A-th  subcarrier  frequency 
using  the  decoding  equations  given  by  [4,  7] 

(A)  =  H*x  (k)rx  [k)  +  Hl2  (A)r3*  (k)  +  H2  l  (A)r2  (A)  +  H2  2  (A)r4*  (A) 


(4.25) 


$2  (A)  =  H*2  (A)fj  (A)-//, ,  (A)r3*  (k)  +  H22  (k)r2  (k)-H2l  (A)r4*  (A) 


(4.26) 


38 


where  Hj  j(k),  i  =  1, 2  j  =  1, 2 ,  are  the  frequency  response  values  of  the  channel  from  the 
j-th  transmit  antenna  to  the  f-th  the  receive  antenna  for  the  A-th  subcarrier,  and 
ri ;  ( A ) ,  i  =  1, 2, 3, 4 ,  are  the  output  of  FFT  blocks  at  the  A-th  subcarrier. 


Figure  19.  Schematic  of  a  2x2  MIMO-OFDM  System  Utilizing  STBC 

C.  SUMMARY 


This  chapter  discussed  the  MIMO  schemes  using  single-carrier  and  OFDM  modu¬ 
lation  over  frequency- selective  channels.  In  the  presence  of  a  frequency  selective  chan¬ 
nel,  STBC  was  used  to  realize  these  MIMO  schemes. 

The  discussion  of  the  single-carrier  and  OFDM  systems  was  based  on  the  simula¬ 
tion  developed  in  Matlab.  While  developing  the  input-output  relations  and  the  decoding 
equations  of  the  systems,  the  transmitted  signals  were  assumed  to  be  in  baseband,  i.e.,  RF 
modulation  was  not  conducted  in  the  simulations.  The  channel  information  given  in  the 
decoding  equations  were  based  on  the  assumption  that  there  was  full  knowledge  of  the 
channel  characteristics  at  the  receiver,  i.e.,  channel  estimation  was  not  performed  in  the 
simulation  process. 
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The  next  chapter  discusses  the  simulation  model  developed  to  study  MIMO-SC 
and  MIMO-  OFDM  schemes  over  the  frequency  selective  SUI  channels  and  presents  the 
simulation  results. 
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V.  SIMULATION  RESULTS 


The  previous  chapters  covered  the  principles  of  the  Alamouti  scheme  over  flat 
channels,  STBC  systems  with  single-carrier  and  OFDM  modulation.  This  chapter  pre¬ 
sents  performance  comparisons  of  the  MIMO  systems  developed  in  Matlab.  In  the  simu¬ 
lation  process,  the  goal  was  to  reach  a  bit  error  rate  (BER)  of  at  least  10  5  to  evaluate  the 
performances  of  the  systems.  Therefore  for  each  Eb/Na  ,  at  least  1,000,000  information 
bits  were  simulated  for  each  system.  One  run  of  sending  about  1,000,000  bits  for  all 
Eb/N0  values  is  called  a  Monte  Carlo  run.  The  systems  were  simulated  for  a  specified 

number  of  Monte  Carlo  runs  and  the  results  were  averaged  such  that  the  BER  curves 
were  smooth  enough  for  evaluation. 

A.  PERFORMANCE  OF  SYSTEMS  USING  ALAMOUTI  SCHEMES 

Initially  a  Matlab  code  was  generated  including  the  SISO  (lxl),  the  SIMO  (1x2), 
the  MISO  (2x1)  and  the  MIMO  (2x2)  systems  with  single-carrier  modulation.  The 
Alamouti  coding  scheme  discussed  in  Chapter  II  was  performed  over  consecutive  symbol 
pairs.  The  developed  systems  were  simulated  over  Rayleigh  fading,  one  tap,  zero  mean, 
circularly  symmetric  and  complex  valued  channels  Hw  (see  properties  of  Hw  in  (3.1)). 
For  the  MISO  and  the  MIMO  cases,  it  was  assumed  that  the  channel  characteristics  re¬ 
mained  constant  over  consecutively  transmitted  symbol  periods.  Each  channel  was  gen¬ 
erated  independently,  i.e.,  it  was  assumed  that  the  channels  in  all  multiple  antenna  cases 
were  uncorrelated.  For  this  simulation,  channel  encoding  was  not  employed.  It  was  as¬ 
sumed  that  the  total  transmitted  power  from  two  antennas  for  the  MISO  and  the  MIMO 
systems  was  the  same  as  the  transmitted  power  from  the  single  antennas  for  the  SISO  and 
SIMO  systems.  Additionally,  we  assumed  that  there  was  full  knowledge  of  the  channel  at 
the  receiver.  After  computing  the  received  signal(s)  at  the  receiver(s),  estimates  of  the 
transmitted  symbols  were  computed  using  the  relations  given  in  (2.13),  (2.17)  and  (2.20). 
Table  7  lists  the  parameters  of  the  simulation  run  using  the  QPSK  modulation.  By  using 
the  QPSK  modulation,  1,000,000  information  bits  were  mapped  into  500,000  symbols. 
The  BER  curves  were  smooth  enough  for  evaluation  with  10  Monte  Carlo  runs. 
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#  of  transmitted  information  bits  for  each  Eb  /  N0 

1,000,000 

#  of  symbols  transmitted  for  each  Eb  /  Nt) 

500,000 

#  of  Monte  Carlo  runs 

10 

Table  7.  Parameters  for  Alamouti  and  MRC  Single-Carrier  Simulations 

Figure  20  shows  the  BER  performance  of  SISO  (lxl),  SIMO  (1x2,  MRC),  MISO 
(2x1,  Alamouti  Scheme)  and  MIMO  (2x2,  Alamouti  Scheme)  systems  over  Rayleigh 
fading  channels.  The  performance  of  the  MISO  (2x1)  is  3  dB  worse  than  the  SIMO 
(1x2),  i.e.,  the  performance  gap  between  the  two  schemes  is  3  dB  at  each  Eb/N0  value. 

This  3  dB  gap  has  occurred  because  the  transmitted  power  for  the  MISO  system  is  kept  at 
half  the  power  for  the  SIMO  system.  If  each  antenna  of  the  MISO  system  had  transmitted 
at  the  same  power  level  as  the  single  transmit  antenna  of  the  SIMO  system,  the  BER 
curves  of  these  two  schemes  would  be  identical.  As  a  result,  we  have  shown  that  a  SIMO 
system  utilizing  the  MRC  scheme  is  identical  to  a  MISO  system  employing  the  Alamouti 
scheme  [4]. 

On  the  other  hand,  both  systems  have  a  significant  improvement  over  the  SISO 
system.  Increasing  the  number  of  antennas  either  at  the  transmitter  or  the  receiver  gave 
spatial  diversity  gain.  Even  at  a  BER  of  2x10  3 ,  the  MISO  system  performs  approxi¬ 
mately  10  dB  better  than  the  SISO  system.  By  observing  the  slope  of  the  curves,  we  can 
remark  that  for  higher  Eb/Na ,  the  gap  will  increase.  Adding  one  more  transmit  antenna 

to  the  SIMO  system  gave  us  more  gain.  Although  the  transmitted  power  was  halved,  the 
MIMO  system  has  about  10  dB  better  performance  than  the  SIMO  system  at  a  BER  of 
10“5 .  As  Eb/Na  is  increased,  the  performance  gap  increased  as  well. 
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Figure  20.  The  BER  Performance  Comparison  of  the  Alamouti  Schemes  to  the  MRC 
Scheme  over  Rayleigh  Fading  Channels  (The  simulation  parameters  are 

given  in  Table  7) 

B.  PERFORMANCE  OF  MIMO-SC  SYSTEMS  OYER  SUI  CHANNELS 

We  discussed  the  SUI  channels  in  Chapter  III  and  the  STBC-SC  schemes  in 
Chapter  IV.  The  parameters  in  STBC-SC  simulations  were  chosen  to  be  comparable  to 
the  OFDM  parameters  of  the  IEEE  802.16a  standard  since  we  wanted  to  compare  the  per¬ 
formances  of  the  single-carrier  schemes  to  the  OFDM  schemes.  To  simulate  the  STBC 
schemes,  a  Matlab  code  was  generated.  The  generation  of  the  SUI  channels  was  achieved 
with  the  modified  version  of  the  Matlab  code  presented  in  [11].  The  output  of  the  channel 
generation  code  was  three  complex  numbers  that  represent  the  three-tap  SUI  channels.  A 
sampling  operation  was  carried  out  to  implement  the  generated  SUI  channels.  The  sam¬ 
pling  frequency  fs  =  22.857  MHz  was  chosen  based  on  the  IEEE  802.16a  standard 

specification.  The  derivation  of  fs  will  be  discussed  later  in  the  OFDM  simulation  sec- 
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tion.  In  STBC,  space  coding  was  performed  over  two  consecutive  blocks  of  symbols  as 
discussed  in  Chapter  III.  To  make  a  fair  comparison  of  the  single-carrier  schemes  to  the 
OFDM  systems,  the  block  size  was  chosen  as  192  since  the  number  of  information  sym¬ 
bols  transmitted  in  an  OFDM  symbol  is  192. 

After  setting  the  sampling  frequency  and  the  symbol  block  size,  discrete-time 
domain  representations  of  the  SUI  channels  were  investigated.  Figure  21  shows  the  three- 
tap  SUI  channels  with  sampling  frequency  fs  =  22.857  MHz  .  This  figure  highlights  the 
different  delay  spreads  encountered  by  different  SUI  channels.  For  example,  the  SUI-3 
channel  model  has  a  second  tap  at  t  =  0.4  //s  and  a  third  tap  t  =  0.9  jus  .  All  other  samples 
were  set  to  zero  in  the  sequence  of  the  sampled  channel.  In  the  implementation,  the  delay 
spread  parameters  matched  the  SUI  channel  parameters  shown  in  Appendix  A.  Notice 
that  SUI-5  and  SUI-6  have  delay  spreads  up  to  10  jus  and  20  jus  ,  respectively. 


time(ji  s) 


time(jj,  s) 


time(|i  s) 


time(fx  s) 


time(|i  s) 


Figure  21.  SUI  Channel  Impulse  Responses  on  a  Sampling  Grid  with  a  Sampling  Fre¬ 
quency  fs  =  22.857  MHz 
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The  next  step  was  to  simulate  the  SISO  system  over  different  SUI  channels.  The 
generation  of  the  192-symbol  block  is  illustrated  in  Figure  22.  The  implementation  util¬ 
izes  a  rate  Vi  convolutional  encoder.  Consequently,  192  information  bits  were  mapped 
into  192  QPSK  symbols.  The  simulation  parameters  for  the  system  using  the  QPSK 
modulation  are  listed  in  Table  8.  In  line  with  the  generation  process  shown  in  Figure  22, 
1,536,000  information  bits  were  mapped  into  4000  symbol  blocks  of  length  192.  In  the 
simulation,  it  was  assumed  that  the  channel  characteristics  remained  constant  over  the 
transmission  period  of  a  192-symbol  block. 


Figure  22.  Schematic  of  Generation  of  192-symbol  Blocks 


#  of  transmitted  information  bits  for  each  Eb  / Na 

1,536,000 

#  of  symbol  blocks  transmitted  for  each  Eh  /  Na 

4000 

Convolutional  Encoder  generator  polynomials 

171,133 

#  of  Monte  Carlo  runs 

30 

Table  8.  Parameters  for  SISO  System  Simulations  over  SUI  Channels 


The  received  signals  and  the  estimates  of  the  transmitted  symbols  were  computed 
using  the  relations  given  in  (4.3)  and  (4.4),  respectively.  The  BER  performance  plots  of 
the  SISO  system  utilizing  QPSK  over  AWGN  and  SUI-1  to  SUI-6  channels  are  shown  in 
Figure  23.  Not  surprisingly  the  system  has  the  best  performance  in  AWGN.  As  the  sys¬ 
tem  experiences  higher  SUI  channels,  the  BER  performance  gets  worse.  The  SUI-3,  4,  5 
and  6  channels  cause  considerable  performance  degradation.  On  the  other  hand,  the  SUI- 
1  and  the  SUI-2  channels  have  less  severe  effects  on  the  system.  In  the  remaining  simula¬ 
tions,  SUI-2  channel  is  used  to  study  the  performance  of  the  MIMO  systems.  SUI-4,  5 
and  6  are  not  further  considered  here  because  of  their  larger  delay  spreads  ( 4  -  20  fis ). 
SUI-2  is  a  reasonable  channel  model  for  an  overall  performance  study  since  it  accounts 
for  multipath  and  yet  has  a  stronger  direct  signal  compared  to  SUI-3. 
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Figure  23.  The  BER  Performance  of  the  SISO-SC  System  over  Different  SUI  Channels 

(The  simulation  parameters  are  given  in  Table  8) 

After  selecting  the  SUI-2  channel  as  the  simulation  environment,  all  STBC  sys¬ 
tems  with  single-carrier  were  simulated  over  the  SUI-2  channel.  For  the  MISO  and  the 
MIMO  systems,  the  192-symbol  blocks  were  generated  as  illustrated  in  Figure  24.  By 
using  a  rate  Vi  convolutional  encoder,  384  information  bits  are  mapped  into  384  QPSK 
symbols,  which  are  then  divided  into  two  symbol  blocks,  one  for  each  antenna.  The  simu¬ 
lation  parameters  of  the  single-carrier  schemes  are  listed  in  Table  9.  As  investigating  the 
performance  of  STBC-SC  schemes  was  one  of  the  major  goals  of  this  work,  the  number 
of  information  bits  and  the  number  of  Monte  Carlo  runs  were  chosen  to  be  higher  than 
those  of  the  previous  simulations.  By  choosing  a  higher  number  of  bits,  we  aimed  to 
measure  the  bit  error  performance  at  lower  rates. 
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192  symbols  to  1st 
antenna 


192  symbols  to  2nd 
antenna 


Figure  24.  Schematic  of  Generation  of  192-symbol  Blocks  for  each  Transmit  Antenna 


#  of  transmitted  information  bits  for  each  Eb  /  Na 

3,840,000 

#  of  symbol  blocks  transmitted  for  each  Eb  /  No 

10,000 

Convolutional  Encoder  generator  polynomials 

171,133 

#  of  Monte  Carlo  runs 

50 

Table  9.  Parameters  for  Simulations  of  the  Single-Carrier  Systems  over  SUI-2 

Channel 


In  the  simulations,  it  was  assumed  that  we  had  the  perfect  channel  knowledge  at 
the  receiver.  We  also  assumed  that  the  channel  characteristics  remained  constant  over 
two  transmission  burst  periods  for  the  MISO  and  MIMO  systems.  The  total  transmitted 
power  is  kept  at  the  same  level  for  all  schemes,  i.e.,  for  the  MISO  and  the  MIMO  cases 
the  power  for  each  antenna  is  halved  as  compared  to  the  SIMO  and  SISO  cases. 

The  equations  used  to  compute  the  received  signals  and  estimates  in  the  simula¬ 
tion  are  (4.3)-(4.4),  (4.6)-(4.7),  (4. 1 1)-(4. 12)  and  (4.25)-(4.26)  for  SISO-SC  (lxl), 
SIMO-SC  (1x2),  MISO-SC  (2x1)  and  MIMO-SC  (2x2)  systems,  respectively. 

The  simulation  results  of  all  single-carrier  schemes  are  shown  in  Figure  25.  The  3 
dB  performance  gap  between  the  MISO  and  the  SIMO  systems  was  observed  again.  The 
MIMO-SC  outperformed  all  other  schemes.  At  a  BER  of  10  4 ,  MIMO-SC  performed 
better  than  the  SIMO-SC  by  about  1.0  dB.  The  performance  gap  is  increased  as  Eh/Na  is 

increased.  The  performance  gap  is  about  2.5  dB  at  a  BER  of  2x10  '' .  The  largest  gap  is 
certainly  between  the  SISO-SC  and  the  MIMO-SC.  At  a  BER  of  6x10  5 ,  the  MIMO-SC 
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outperforms  the  SISO-SC  by  approximately  11  dB.  This  result  showed  that  adding  more 
antennas  at  both  ends  of  the  communication  link  gave  us  significant  performance  im¬ 
provement. 
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Figure  25.  The  BER  Performance  Comparisons  of  Single-Carrier  Schemes  over  the 
SUI-2  Channel  (The  simulation  parameters  are  given  in  Table  9) 


The  next  simulation  was  conducted  to  observe  the  effects  of  spatial  correlation  on 
the  single-carrier  systems.  The  SISO  and  the  MISO  systems  were  chosen  for  examining 
the  correlation  effects.  The  BER  performance  of  the  MISO  system  with  different  transmit 
correlation  coefficients  pt  is  compared  to  that  of  the  SISO  case.  The  transmit  correlation 

matrix  was  constructed  as  given  in  (3.5).  The  correlated  MISO  channel  was  rebuilt  with 
the  relation  given  in  (3.8).  The  simulation  parameters  of  the  MISO  system  over  the  corre¬ 
lated  channel  are  listed  in  Table  10.  By  using  the  QPSK  modulation,  1,536,000  informa¬ 
tion  bits  were  mapped  into  8000  symbols  blocks,  each  of  length  192.  The  BER  curves 
were  plotted  based  on  30  Monte  Carlo  runs. 
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#  of  transmitted  information  bits  for  each  Eh  / Na 

1,536,000 

#  of  symbol  blocks  transmitted  for  each  Eh  / Na 

8,000 

Convolutional  Encoder  generator  polynomials 

171,133 

#  of  Monte  Carlo  runs 

30 

Table  10.  Parameters  for  Simulations  of  the  MISO  System  over  Correlated  SUI-2 

Channels 


The  BER  performance  of  the  MISO  single-carrier  system  over  correlated  channels 
is  shown  in  Figure  26.  At  a  BER  of  2xl0~4 ,  the  MISO  system  performed  about  4  dB  bet¬ 
ter  than  the  SISO  system  over  uncorrelated  channels.  As  the  correlation  coefficient  in¬ 
creased,  the  BER  got  worse.  For  example,  when  pt  =  0.9  the  performance  gap  at  the 

BER  of  2x10^  decreased  to  2  dB.  As  a  result,  correlation  between  the  channels  caused 
significant  performance  degradation  on  the  systems  with  multiple  antennas.  The  correla¬ 
tion  coefficients  pt  and  pr  mostly  depend  on  the  spatial  parameters,  such  as  distance 

between  the  antennas  and  the  wavelength.  Therefore,  placing  the  antennas  and  selecting 
the  wavelength  play  important  roles  on  the  performance  of  the  systems. 

C.  PERFORMANCE  OF  MIMO-OFDM  SYSTEMS  OYER  SUI  CHANNELS 
The  final  goal  of  this  thesis  was  to  investigate  the  performance  of  MIMO-OFDM 
systems.  The  OFDM  parameters  were  chosen  from  the  IEEE  802.16a  standard.  The  pa¬ 
rameters  in  the  IEEE  802.16a  standard  were  specified  for  various  bandwidths.  For  this 
thesis,  the  bandwidth  (W)  was  chosen  as  20  MHz  for  a  license-exempt  band.  The  sam¬ 
pling  frequency  fs  was  derived  from  the  ratio  [10] 


W  7  ' 


(5.1) 


Substituting  W  =  20  MHz  into  (5.1),  the  sampling  frequency  was  determined  to  be 


fs  =  22.857  MHz  . 


(5.2) 
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Figure  26.  The  BER  Performance  of  the  MISO  Single-Carrier  System  over  Correlated 
Channels  (The  simulation  parameters  are  given  in  Table  10) 


The  sampling  frequency  is  the  most  important  parameter  in  the  simulation,  be¬ 
cause  the  sampling  of  the  impulse  response  of  the  three-tap  SUI  channel,  subcarrier  spac¬ 
ing,  useful  time,  guard  interval  and  the  symbol  duration  are  all  based  upon  it.  The  carrier 
spacing  was  derived  from  the  relation 


A/  = 


/, 

N 

1  v  FFT 


(5.3) 


where  NFFT  is  the  FFT  size.  Substituting  NFFT  =  256  and  fs  =  22.857  MHz  ,  the  carrier 


spacing  was  calculated  to  be 

A/  =  89.3  kHz  . 

The  useful  time  Tb  was  computed  with  the  relation 

T,=  — =  11.2  {is . 
Af 


(5.4) 


(5.5) 
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The  total  OFDM  period  is  given  to  be 

Ts=Tg+Tb  (5.6) 

where  T  is  the  guard  interval  in  an  OFDM  symbol.  The  guard  interval  is  determined 
from  the  optional  guard  ratios  T/Tb  .  Table  11  lists  the  OFDM  symbol  durations  for 
various  guard  ratios,  which  range  for  1/32  to  1/4.  By  choosing  the  smallest  guard  ratio 
Tg /Tb  =  1/32 ,  the  total  OFDM  symbol  duration  T  is  pretty  close  to  the  useful  time 

Tb  =11.2  /us  .  As  the  guard  ratio  T  /Tb  increases,  the  length  of  CP  and  OFDM  symbol 

also  increases.  The  guard  time  (i.e.,  the  CP  length)  is  chosen  to  be  larger  than  the  ex¬ 
pected  delay  spread  of  the  channel  so  that  multipath  components  from  one  symbol  would 
not  interfere  with  those  from  the  next  symbol.  In  choosing  the  guard  time,  we  trade-off 
between  robustness  against  multipath  and  bandwidth  efficiency.  To  mitigate  the  effects  of 
multipath,  a  large  guard  time  can  be  chosen,  which,  on  the  other  hand,  it  requires  wider 
bandwidth. 


Guard  Ratio 

CP  length  with  respect  to 

OFDM  symbol 

T./T> 

the  FFT  size  NFFT  =256 

duration  Ts 

1/32 

8 

11.55  /us 

1/16 

16 

11.9  jus 

1/8 

32 

12.6  jus 

1/4 

64 

14.0  jus 

Table  11.  OFDM  Symbol  Durations  for  Various  Guard  Ratios 


After  selecting  the  OFDM  channelization  parameters,  a  SISO-OFDM  system  was 
simulated  and  the  performances  were  compared  to  those  of  the  SISO-SC  system.  The 
SISO-OFDM  system  was  simulated  without  CP  and  with  the  largest  guard  ratio 
Tg  jTb  =  1/4 ,  i.e.,  from  Table  1 1  the  CP  size  is  64.  The  simulation  parameters  are  listed  in 


51 


Table  12.  By  using  the  QPSK  modulation,  1,536,000  information  bits  were  mapped  into 
8,000  symbols  blocks  of  length  192  and  8,000  OFDM  symbols. 


#  of  transmitted  information  bits  for  each  Eh  /  Na 

1,536,000 

#  of  symbol  blocks  transmitted  for  each  Eh  /  Na 

8,000 

#  of  OFDM  symbols  transmitted  for  each  Eb/Na 

8,000 

Convolutional  Encoder  generator  polynomials 

171,133 

#  of  Monte  Carlo  runs 

30 

Table  12. 


Parameters  for  Simulations  of  the  SISO-SC  and  SISO-OFDM  Systems 
over  the  SUI-2  Channel  Model 


The  BER  comparison  of  SISO-OFDM  to  the  SISO-SC  is  shown  in  Figure  27.  The 
SISO-OFDM  system  without  CP  has  about  1  dB  better  performance  at  the  BER  of  10  3 
when  compared  to  the  SISO-SC.  For  lower  bit  error  rates,  the  performance  of  the  SISO- 
OFDM  gets  much  better.  At  a  BER  of  10  5 ,  the  performance  difference  is  about  6  dB  al¬ 
though  no  CP  was  added  at  the  transmitter.  This  result  is  due  to  the  SISO-OFDM  sys¬ 
tem’s  robust  performance  over  frequency  selective  channels  compared  to  the  SISO-SC 
scheme.  When  the  guard  ratio  was  chosen  as  1/4 ,  the  BER  curve  moved  to  left  about  1 
db. 

The  last  step  was  to  present  all  multiple  antenna  schemes  in  one  graph.  To  make  a 
fair  comparison,  the  symbol  block  size  of  the  single-carrier  scheme  was  selected  to  be 
192.  The  guard  ratio  Tg  jTh  in  all  the  OFDM  schemes  was  chosen  to  be  1/4 .  For  all  mul¬ 
tiple  transmit  antenna  schemes  (MISO-SC,  MISO-OFDM,  MIMO-SC,  MIMO-OFDM), 
the  transmitted  power  from  each  antenna  is  halved  when  compared  to  the  single  transmit 
antenna  schemes  (SISO-SC,  SISO-OFDM,  SIMO-SC,  SIMO-OFDM).  The  channel  char¬ 
acteristics  were  assumed  to  be  constant  over  consecutive  bursts  of  symbol  blocks  or 
OFDM  symbols.  The  multiple  channels  were  generated  independently,  i.e.,  it  was  as¬ 
sumed  that  all  channels  were  uncorrelated. 
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Figure  27.  The  BER  Comparison  of  SISO-OFDM  with  Various  CP  Sizes  to  SISO-SC 
over  the  SUI-2  Channel  (The  simulation  parameters  are  given  in  Table  12) 

The  equations  used  to  compute  the  received  signals  and  estimates  in  the  simula¬ 
tion  are  (4.18)-(4.19),  (4.20)-(4.21),  (4.23)-(4.24)  and  (4.25)-(4.26)  for  SISO-OFDM 
(lxl),  SIMO-OFDM  (1x2),  MISO-OFDM  (2x1)  and  MIMO-OFDM  (2x2)  systems,  re¬ 
spectively. 

The  simulation  parameters  for  all  schemes  are  listed  in  Table  13.  The  parameters 
were  chosen  to  be  of  the  same  as  in  the  single-carrier  simulation  (see  Table  9).  By  choos¬ 
ing  the  QPSK  simulation,  3,840,000  information  bits  were  mapped  into  10,000  OFDM 
symbols  and  192-symbol  blocks. 

The  BER  plots  for  all  schemes  are  shown  in  Figure  28.  The  dashed  lines  are  for 
the  SC  schemes  and  the  solid  lines  are  for  the  OFDM  schemes.  Obviously  the  MIMO- 
OFDM  system  outperforms  all  other  schemes  by  reaching  a  BER  of  10  5  with 
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Eh/Na  =  7.5  dB  .  When  the  MIMO-OFDM  is  compared  to  the  worst  performing  system, 

SISO-SC,  the  performance  difference  is  about  12  dB  at  a  BER  of  10  4 .  The  3-dB  gap  is 
again  observed  between  the  SIMO-OFDM  and  the  MISO-OFDM  systems.  Although 
there  are  totally  three  antennas  used  in  the  SIMO-OFDM  system,  the  SIMO-OFDM  per¬ 
forms  better  than  MIMO-SC  (totally  four  antennas)  until  the  BER  is  about  10~5 .  The  per¬ 
formance  of  the  two  schemes  is  fairly  close  after  this  BER  crossover  point. 


#  of  transmitted  information  bits  for  each  Eh  /  Na 

3,840,000 

#  of  symbol  blocks  transmitted  for  each  Eb  /  No 

10,000 

#  of  OFDM  symbols  transmitted  for  each  Eb/Na 

10,000 

Convolutional  Encoder  generator  polynomials 

171,133 

#  of  Monte  Carlo  runs 

50 

Table  13.  Parameters  for  Simulations  of  the  all  Single-Carrier  and  OFDM  Systems 

over  the  SUI-2  Channel 

D.  SUMMARY 

This  chapter  has  presented  the  simulation  results  of  several  Matlab-based  com¬ 
munication  systems.  Simulations  were  performed  over  three  main  steps.  First,  the 
Alamouti-scheme -based  MIMO  systems  were  simulated  over  Rayleigh  fading,  single  tap, 
zero  mean  and  circularly  symmetric  channels.  Next,  the  simulations  of  all  single-carrier 
systems  over  frequency-selective  SUI  channels  were  conducted.  Fastly,  all  OFDM  sys¬ 
tems  were  simulated.  As  expected,  overall  the  systems  with  more  antennas  performed 
better  than  the  systems  with  fewer  antennas.  Another  expected  result  was  that  the  OFDM 
systems  performed  better  than  the  systems  utilizing  single-carrier  modulation. 
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BER 


Figure  28. 


The  BER  Comparison  of  all  Schemes  over  the  SUI-2  Channel  (The  simula¬ 
tion  parameters  are  given  in  Table  13) 
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VI.  CONCLUSION 


The  goal  of  this  thesis  was  to  compare  the  performance  of  MIMO  systems  to  the 
conventional  SISO  systems.  This  goal  was  accomplished  by  investigating  the  STBC 
technique  and  the  corresponding  decoding  algorithms.  The  investigated  space  time  cod¬ 
ing  techniques  were  combined  with  OFDM  to  develop  the  MIMO-OFDM  systems  and 
the  BER  performance  of  these  systems  were  obtained  through  a  comprehensive  simula¬ 
tion  study. 

A.  SUMMARY  OF  THE  WORK  DONE 

A  comprehensive  background  of  the  MIMO  systems  was  presented,  including  the 
input-output  relations,  space-coding  techniques  over  flat  and  frequency  selective  channels 
and  the  decoding  algorithms.  The  presented  schemes  were  investigated  for  both  the  sin¬ 
gle-carrier  and  the  OFDM  modulation  cases.  For  the  simulation  environment,  the  fre¬ 
quency-selective  SUI  channels  were  utilized. 

Matlab-based  MIMO-SC  and  MIMO-OFDM  communication  systems  were  de¬ 
veloped.  The  simulations  were  performed  for  both  uncorrelated  and  correlated  MIMO 
channels  based  on  the  SUI  models.  Simulation  results  were  presented  in  the  form  of  the 
BER  performance  curves  and  compared  among  the  various  MIMO  systems  studied  in  this 
work. 

B.  SIGNIFICANT  RESULTS  AND  CONCLUSIONS 

The  most  significant  result  of  this  study  was  the  observation  that  the  system  per¬ 
formance  improved  as  the  number  of  antennas  at  either  side  of  the  communication  system 
was  increased.  As  more  antennas  were  added  to  the  conventional  single  antenna  systems, 
higher  diversity  gains  improved  the  performance.  Obviously,  the  system  with  the  largest 
number  of  antennas,  i.e.,  the  2x2  system  in  this  work,  outperformed  all  other  systems 
with  fewer  antennas  by  exploiting  the  spatial  diversity  gain. 

The  3-dB  penalty  for  the  MISO  systems  utilizing  STBC  when  compared  to  the 
SIMO  systems  utilizing  MRC  was  observed  in  all  single-carrier  and  OFDM  schemes. 

This  penalty  has  occurred  because  the  transmitted  power  for  each  antenna  was  halved  for 
the  MISO  systems  when  compared  to  the  transmitted  power  for  the  SIMO  systems.  On 
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the  other  hand,  if  the  total  transmitted  power  is  doubled  for  the  MISO  system,  then  its 
performance  was  found  to  be  identical  to  the  SIMO  system.  In  general,  building  multiple 
antennas  at  the  base  stations  (i.e.,  at  the  transmitter)  is  more  feasible  than  building  multi¬ 
ple  antennas  at  the  mobile  stations  (i.e.,  at  the  receiver)  [4].  Moreover,  to  achieve  the 
same  performance,  having  multiple  antennas  at  the  base  station  is  much  more  profitable 
than  having  multiple  antennas  at  all  mobile  stations. 

The  other  significant  result  was  that  the  combination  of  OFDM  and  MIMO 
yielded  significant  performance  improvement  over  frequency- selective  multipath  chan¬ 
nels.  Increasing  the  number  of  antennas  at  both/either  the  receiver  and/or  the  transmitter 
added  more  gain  to  the  OFDM  systems.  The  use  of  multiple  antennas  added  spatial  diver¬ 
sity  gain  to  the  conventional  single-antenna  OFDM  systems.  Therefore,  MIMO-OFDM 
systems  appear  to  be  attractive  for  future  wireless  communication  systems. 

C.  SUGGESTIONS  FOR  FUTURE  WORK 

There  are  four  areas  identified  for  future  work.  First,  the  channel  information  can 
be  extracted  from  the  pilot  carriers  at  the  receiver.  In  all  simulations,  we  assumed  that 
there  was  perfect  channel  knowledge  at  the  receiver.  Especially  for  the  OFDM  simula¬ 
tions,  having  known  the  channel  (i.e.,  all  three  taps  of  the  SUI  channels)  helped  us  to 
compute  the  exact  frequency  response  of  the  channel  at  all  256  subcarrier  frequencies. 
Therefore,  in  the  decoding  process  for  all  schemes,  perfect  channel  knowledge  boosted 
the  performance.  In  practice,  the  channel  information  must  be  estimated  by  using  the  pi¬ 
lot  subcarriers.  Accordingly,  the  frequency  response  of  the  channel  is  known  at  the  re¬ 
ceiver  with  some  inaccuracies.  These  inaccuracies  may  lead  to  performance  degradation. 
In  this  work,  while  computing  the  IFFT  in  the  Matlab  code,  the  pilot  symbols  were  set  to 
null.  These  null  spaces  may  be  filled  with  pilot  symbols  and  the  performance  of  the 
MIMO  systems  may  be  investigated  with  the  channel  information  extracted  from  the  pilot 
subcarriers. 

A  more  extensive  channel  encoder  may  be  used.  The  specifications  for  the  OFDM 
systems  were  chosen  from  the  IEEE  802.16a  standard.  Although  the  concatenated  FEC 
encoders  or  block  turbo  encoders  were  introduced  in  the  standard,  only  a  rate  Vi  convolu- 
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tional  encoder  was  used  in  all  simulations  in  this  work.  Therefore,  the  Matlab  code  may 
be  enhanced  with  a  concatenated  channel  encoder. 

In  the  correlated  MIMO  channels,  the  correlation  coefficients  pt  and  pr  were 

chosen  arbitrarily  in  this  thesis.  In  the  literature,  several  authors  investigated  the  corre¬ 
lated  MIMO  channels  [14,  19-22].  New  Matlab  code  may  be  generated  to  compute  the 
correlation  coefficients  with  respect  to  the  antenna  spacing,  wavelength,  the  angle  of  arri¬ 
val  and  the  angle  of  departure  of  the  electromagnetic  waves. 

Lastly,  the  number  of  antennas,  either  at  the  receiver  or  at  the  transmitter,  may  be 
increased.  The  space  coding  techniques  for  higher  number  of  antennas  based  on  the 
Alamouti  scheme  were  investigated  by  [23,  24].  With  higher  number  of  antennas,  the  de¬ 
coding  equations  will  be  much  more  complicated  but  performance  improvements  may  be 
achieved  due  to  possible  spatial  diversity  gain. 
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APPENDIX  A.  SUI  CHANNEL  PARAMETERS 


This  appendix  lists  the  SUI  channel  parameters.  Each  channel  model  has  different 
parameters  for  the  omni  and  30°  directional  antenna  systems.  The  2nd  tap  power  of  the 
30  directional  channel  is  attenuated  by  6  dB  and  the  3  power  is  attenuated  by  an  addi¬ 
tional  12  dB  compared  to  the  omnidirectional  cases,  /('-factors  are  defined  for  90%  and 
75%  cell  coverage.  Only  the  90%  coverage  parameters  are  shown  here.  More  detailed 
channel  parameters  are  available  in  [1 1]. 


SUI-1  Channel 

Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

0.4 

0.9 

ps 

Power  (omni  ant.) 

0 

-15 

-20 

dB 

Power  (30  ant.) 

0 

-21 

-32 

dB 

/Cfactor  (omni  ant.) 

4 

0 

0 

dB 

K  -factor  (30  ant.) 

16 

0 

0 

dB 

Doppler 

0.4 

0.3 

0.5 

Hz 

SUI-2  Channel 

Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

0.4 

1.1 

ps 

Power  (omni  ant.) 

0 

-12 

-15 

dB 

Power  (30  ant.) 

0 

-18 

-27 

dB 

K  -factor  (omni  ant.) 

2 

0 

0 

dB 

K  -factor  (30  ant.) 

8 

0 

0 

dB 

Doppler 

0.2 

0.15 

0.25 

Hz 

SUI-3  Channel 

Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

0.4 

0.9 

ps 

Power  (omni  ant.) 

0 

-5 

-10 

dB 

Power  (30  ant.) 

0 

-11 

-22 

dB 

K  -factor  (omni  ant.) 

1 

0 

0 

dB 

K  -factor  (30  ant.) 

3 

0 

0 

dB 

Doppler 

0.4 

0.3 

0.5 

Hz 
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SUI-4  Channel 


Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

1.5 

4 

ps 

Power  (omni  ant.) 

0 

-4 

-8 

dB 

Power  (30  ant.) 

0 

-10 

-20 

dB 

K  -factor  (omni  ant.) 

0 

0 

0 

dB 

K  -factor  (30  ant.) 

1 

0 

0 

dB 

Doppler 

0.2 

0.15 

0.25 

Hz 

SUI-5  Channel 

Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

4 

10 

ps 

Power  (omni  ant.) 

0 

-5 

-10 

dB 

Power  (30  ant.) 

0 

-11 

-22 

dB 

K  -factor  (omni  ant.) 

0 

0 

0 

dB 

K  -factor  (30  ant.) 

0 

0 

0 

dB 

Doppler 

2 

1.5 

2.5 

Hz 

SUI-6  Channel 

Tap-1 

Tap- 2 

Tap- 3 

Units 

Delay 

0 

14 

20 

ps 

Power  (omni  ant.) 

0 

-10 

-14 

dB 

Power  (30  ant.) 

0 

-16 

-26 

dB 

K  -factor  (omni  ant.) 

0 

0 

0 

dB 

K  -factor  (30  ant.) 

0 

0 

0 

dB 

Doppler 

0.4 

0.3 

0.5 

Hz 

Table  14.  Parameters  of  the  SUI  Channel  Models 
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APPENDIX  B.  MATLAB  CODE  EXPLANATION 


This  appendix  briefly  explains  the  Matlab  code,  the  function  of  each  m-file,  and 
the  interrelationships. 

To  simulate  the  discussed  schemes,  two  outer  functions  were  generated.  The  func¬ 
tion  alamouti_schemes_over_rayleigh.m  was  used  for  simulating  the  single-carrier 
Alamouti  schemes  over  Rayleigh  fading  channels.  For  simulating  the  single-carrier  and 
the  OFDM  schemes  over  the  SUI  channels,  the  function  a ll_ofdm_and_single_carrier 
_schemes. m  was  used. 

A.  SIMULATION  PARAMETERS 

Each  function  starts  by  reading  the  simulation  parameters.  The  simulation  pa¬ 
rameters  are  divided  into  four  sections:  main  parameters,  channel  parameters,  noise  pa¬ 
rameters,  and  channel  encoder  parameters. 

The  variables  in  each  parameter  section  are  listed  in  Table  15.  The  total  number 
of  information  bits  to  be  transmitted  for  the  given  Eh  /  Na  value  can  be  chosen  such  that 
the  resulting  BER  rates  are  satisfactory  for  evaluation.  In  the  simulations,  it  was  chosen 
to  be  at  least  1,000,000  to  reach  bit  error  rates  on  the  order  of  10  5 .  For  the  OFDM  and 
single-carrier  simulations  over  SUI  channels,  it  was  chosen  to  be  a  multiple  of 
192xlog2  ( M )  to  match  the  192-symbol  block  size. 

The  number  of  Monte  Carlo  runs  were  controlled  by  the  variable  ‘MonteCarlo’ . 
Before  choosing  ‘MonteCarlo’,  it  is  recommended  that  the  Profiler  function  of  Matlab  be 
used.  By  setting  up  the  ‘MonteCarlo’  as  1,  the  total  time  of  the  desired  simulation  can  be 
measured.  Then,  the  expected  time  of  the  simulation  can  be  estimated  easily,  for  larger 
number  of  Monte  Carlo  runs. 
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Main  Parameters 

N 

FFT  size 

M 

m-ary  PSK  constellation 

N 

total  number  of  information  bits  to  be  transmitted  for  the  given 

Eb/N0  value 

Tg_rate 

CP  rate  with  respect  to  the  FFT  size 

MonteCarlo 

number  of  Monte  Carlo  runs  to  be  simulated 

EbNo_db 

Eh/Na  interval  to  be  simulated  (in  dB) 

Channel  Parameters 

SUI_index 

SUI  model  index,  from  1-6 

direct 

antenna  directivity,  omni  antenna  =  0,  30  deg  directional  =  1 

rho_tx 

pt  transmit  correlation  coefficient 

rho_rx 

pr  receive  correlation  coefficient 

fs 

sampling  frequency 

Noise  Parameters 

sig 

standard  deviation  of  the  noise  sequence 

No 

linear  noise  power 

No_db 

noise  power  in  dB 

Convolutional  Encoder  Parameters 

constlen 

constraint  length 

code gen 

generator  polynomials  in  octal  form 

trel 

trellis  description  of  the  encoder 

Table  15.  The  Variables  Representing  the  Simulation  Parameters 
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B.  CONDUCTED  ITERATIONS 

After  the  system  is  configured  with  these  parameters,  three  main  iterations  start.  A 
conceptual  diagram  of  the  iterations  and  the  steps  taken  in  each  iteration  are  shown  in 
Figure  29.  The  main  iteration  is  the  Monte  Carlo  iteration  controlled  by  the  variable 
‘MonteCarlo’.  The  next  iteration  is  the  Eh/N()  iteration  controlled  by  the  variable 

‘EbNo_db\  The  most  inner  loop  is  the  bit  iteration  controlled  by  the  variable  ‘n’.  From 
Figure  29,  all  system  steps  are  conducted  in  the  bit  iteration.  In  other  loops,  such  as 
Monte  Carlo  and  Eh/Na  iterations,  the  number  of  bit  errors  are  summed  up  and  saved  for 

computing  the  bit  error  rate  subsequently.  The  numbers  of  bit  errors  are  saved  until  the 
Monte  Carlo  loop  ends.  To  compute  the  BER,  the  total  number  of  errors  is  simply  di¬ 
vided  by  the  iterated  number  of  bits. 
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•  compute  the  BER  dividing  the  total  number  of  errors  by  the 
simulated  number  of  bits 

Figure  29.  A  Conceptual  Diagram  of  the  Iterations  for  the  Outer  Functions 


C.  EXPLANATIONS  OF  THE  SUB-FUNCTIONS 

While  generating  the  outer  codes,  both  Matlab’s  built-in  and  user  constructed 
functions  were  used.  Here,  we  explain  the  important  functions  used  in  the  code. 

randint.m  is  the  Matlab’s  built-in  function  generating  random  bits  "0"  or  "1"  with 
equal  probability. 


66 


convene. m  is  the  Matlab’s  built-in  function  which  convolutionally  encodes  the 
randomly  generated  bits.  The  configuration  of  the  convolutional  encoder  can  be  changed 
in  the  channel  encoder  section. 

bin_2_mary.m  converts  the  convolutionaly-encoded  bits  to  m- ary  symbols.  The 
resulting  sequence  of  this  function  is  a  vector  with  decimal  numbers  representing  the 
symbols  in  the  given  constellation. 

clmodce.m  is  the  Matlab’s  built-in  function  which  makes  the  m- ary  mapping  op¬ 
eration.  Decimal  numbers,  representing  the  symbols  in  the  constellation,  are  turned  into 
complex  numbers. 

ifft_802_16.m  takes  a  192-symbol  block  sequence  and  creates  an  OFDM  symbol. 
The  function,  organizes  the  information  symbols  as  specified  in  the  IEEE  802.16a  stan¬ 
dard  (see  Figure  15  and  Table  6).  Prior  to  the  IFFT  operation,  the  guard  and  pilot  subcar¬ 
riers  are  set  to  zero.  After  the  IFFT  operation,  the  function  copies  the  desired  CP  length 
elements  and  adds  them  to  the  beginning  of  the  output  of  the  IFFT  block.  The  CP  length 
can  be  controlled  by  the  variable  ‘Tg_rate’  which  is  the  optional  guard  ratio  Tg/Tb  .  The 

guard  ratio  can  be  chosen  as  1/4,1/8,1/16  or  1/32  .  The  resulting  CP  length  will  be  pro¬ 
portional  to  the  FFT  size  256.  The  output  of  the  function  is  the  baseband  representation  of 
an  IEEE  802.16a  OFDM  symbol. 

SUI_model.ni  generates  a  3x1  complex  vector  including  three-tap  values  of  a  SUI 
channel.  All  six  SUI  channel  parameters  are  saved  in  this  function.  The  function  gener¬ 
ates  the  channel  tap  values  using  the  desired  SUI  index  (1  to  6)  and  the  directivity  of  the 
antenna  (omni  or  30°  directional  antenna). 

SUl_model_resample.m  holds  the  delay  spread  values  for  all  six  SUI  channels. 

The  function  takes  a  3  tap  SUI  channel  and  samples  it  with  the  desired  sampling  fre¬ 
quency. 

create _mimo_channel.m  generates  the  channels  for  a  2x2  system.  The  function 
calls  SUI_model.ni  four  times  and  generates  the  uncorrelated  four  paths  of  the  2x2 
MIMO  channel  matrix.  Next,  the  correlated  channel  is  computed  using  the  fading  com- 
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ponents  of  the  four  SUI  channels  and  the  desired  correlation  coefficients  ( pt  and  pr ).  The 
complete  correlated  MIMO  channel  is  computed  with  adding  the  fading  component  to  the 
fixed  component.  The  diagram  shown  in  Figure  9  gives  the  steps  taken  in  the  function. 
Lastly,  the  function  SU I _model_re sample. m  is  called  to  sample  the  four  MIMO  paths 
with  the  desired  sampling  frequency. 

create _miso_channel.m  generates  the  channels  for  a  2x1  system.  The  steps  taken 
in  generating  the  correlated  MISO  channel  are  similar  to  those  for  the  function  cre¬ 
ate  _mimo_channel.m. 

create _simo_channel.m  generates  the  channels  for  a  1x2  system.  The  steps  taken 
in  generating  the  correlated  SIMO  channel  are  similar  to  those  for  the  function  cre¬ 
ate  _mimo_channel.m. 

add_noise.m  generates  a  sequence  of  complex  random  numbers  chosen  from  a 
normal  distribution  with  zero  mean  and  desired  noise  variance.  The  dimension  of  the 
generated  sequence  is  the  same  as  the  sequence  that  noise  will  be  added. 

j ft_802_16.nl  takes  the  received  OFDM  symbol,  removes  the  CP,  takes  the  FFT 
and  organize  the  information  carriers.  The  output  is  a  sequence  with  dimensions  1x192. 

freq_res_info_carriers.m  computes  the  frequency  response  of  the  corresponding 
192  information  subcarriers.  The  information  subcarrier  allocation  is  identical  to  the 
functions  ifft_802_16.m  and  fft_802_  1 6. m 

non_causal_filter.m  is  used  in  decoding  equations  when  a  non-causal  filter  opera¬ 
tion  is  needed. 

ddemodce.m  is  the  Matlab’s  built-in  function  which  makes  the  m- ary  demapping 
operation.  The  function  takes  the  computed  estimated  complex  sequences  and  demaps 
them  to  decimal  numbers,  representing  the  symbols  in  the  constellation. 

mary_2_bin.m  converts  the  decimal  number  to  bits  by  using  the  desired  constella¬ 
tion  size. 

Vitdec.m  is  the  Matlab’s  built-in  function  which  makes  the  Viterbi  decoding  algo¬ 
rithm.  In  the  simulation,  the  function  is  configured  for  hard  decision  decoding. 
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